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ABSTRACT 


The BRRVAR procedure provides estimates of the samp!ing error or sta- 
tistics generated from complex samples. BRRVAR uses the Balanced Repeated 
Replication approach (BRR), which requires that the primary sampling units 
be organized into pairs that reflect the act:uai sampling design, Common 
statistics are generated, including sums, means, standard deviations, co- 
variances, correlations, and regression coefficients. For each statistic, 
an overall estimate is computed along with the estimate of the standard 
error of the statistic and the appreximate 45% confidence bounds for the 


statistic. 
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INTRODUCTION 


BRRVAR is a user-created procedure for use with the Statistical Analysis 
System (SAS). It was created for tle National Center for Education Statistics 
to enlarge their capacity to estimate and analyze sampling errors for :tatis~ 
tics generated from educational surveys with complex sampling designs. ‘This 
program drew heavily from an earlier SAS procedure, NASSVAR, desiguued by 
David Morganstern and implemented by Greg Binzer, both of Westat, Inc. The 
most significant enhancements in the current program are the “automatic” 
application of an appropriate design matrix for defining the replicate half- 


samples and the much wider range of statistics which can be requested simply. 


The use of BRRVAR may require considerable preparation of the input data 
file. If the initial sampling design was created specifically for BRR ana- 
lyses, it specified sampling strata from which exactly two units were randomly 
drawn. (If the sampling was multi-stage, e.g. schools were sampled and then 
students were sampled within school, only the first or primary stage is of 
imporance here.) In this case, all that is required is that the file contain 
one or more variables that can by used to identify the different strata and 
the different PSUs within each stratum. The file must be sorted by these 
variables prior to calling BRRVAR. 


In most cases, the actual sample drawn does not match the BRR assumptions 
exactly. It is then necessary to create pairs of primary sampling units and 
sort the tile so that these pairs are adjacent. In some instances, the number 
of primary sampling units may be quite large. This is true, for example, if 
some primary sampling units were sampled with certainty so that the secondary 
units within these primary sampling units have to be treated as if they were 
the primary units. In such cases, it may be desirable to define "pseudo-PSUs" 
which are actualiy clusters of PSUs. If several schools were sampled with 
certainty for example, and then some number of students were sampled from each 
school, it would be desirable to define two randomly equivalent clusters of 
students for each school and then treat these clusters as "pseudo~PSUs". 
Similarly, if a relatively large number of schools were selected from some 
sampling stratum, it would be desirable to organize these schools into two 


randomly equivalent clusters and then use these clusters as “pleudo-PSUs." 
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In addition to specifiyng the variables that identify primary sampling 
units, it is necessary to specify a case weight using the WEIGHT statement. 
It is assumed that such a weight will have been put on the file for data 
from the kind of complex samples for which BRRVAR estimates are required. 

In general, the sampling weight will be the inverse of the overall sampling 
probability for each case (particularly if sums are of interest). If no such 


weight exists, then a dummy weight must be created. 
Sepcification 
The following statements are used with the BRRAVAR procedure: 


PROC BRRVAR options and parmms; 
VAR(IABLES) variable list; 
WEIGHT weight variable; 
STRATID variable list; 


PSULD variable list. 


The PROC BRRVAR statement must include certain required parameters as 
described below. The VARIABLES, WEIGHT, STRATID, and PSUID statements are 


also required. Their function and syntax is described below. 


PROC BRRVAR Statement 


The following options may be provided in the PROC BRRVAR statement: 


® DATA = <SAS dataset>, names the SAS dataset to be used by 
PROC BRRVAR. If no input dataset is specified, 
BRRVAR uses the most recently created SAS dataset. 

@ SUMS = causes sums to be estimated for each variable in 
the VARIABLE list. 

@ MEAN = causes means to be estimated for each variable in 


the VARIABLE list. 
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STD 


cov 


CORR 


REGR 


NOMISS 


DESIGN 


DUMP 


NPSU 


NREP 


Causes standard deviations to be estimated for each 
variable in the variable list. 


causes covariances to be estimated for each varia- 
ble in the variable list. (Note: this option is 
not compatible with the CORR option; if both are 
specified, correlations take precedence and cova- 
riances will not be printed.). 


causes correlations to be estimated for each pair 
of variables in the VARIABLE list. 


causes regression coefficients and a multiple R 

square to be estimated for each regresson equation 
specified. The last NDEP variables (1 if NDEP is 

not specified are treated as criterion variables. 

The remaining variables are treated as predictors. 

A separate regression equation is estimated for each 
criterion, using all of the predictor variables for 
each equation. If CORR has been specified, standar- 
dized regresssion coefficients are computed, otherwise 
raw regression coefficients are computed. 


causes listwise deietion of missing values (i.e., if 
any variable is missing, the entire case is excluded 
from all statistics). If only univariate statistics 
are requested and NOMISS is not specified, missing 

values will be deleted separately for each variable. 
If any of the multivariate statistics are requested 

(COV, CORR, or REGR) NOMISS is automatically forced. 


causes the replicate half-sample design matrix to be 
printed. The design matrix contains one row for each 
PSU aad one column for each half-sample. The entries 
tell whether the PSU associated with the row was 
included in the half-sample indicated by the column. 


requests a printout of the half-sample estimates for 
each statistic. In the current version this dump is 
relatively unformatted, with no labelling of the 
columns or rows. 


<number of PSUs>, specifies the number of PSUs in the 
sample design. This parameter is primarily for docu- 
mentation purposes at present. The actual number of 
PSUs is computed and compared to the expected number 
specified here. In case of discrepancies, the acutal 
Number is used. 


<number of replicates>, specifies the number of half- 
sample replicates to be used. If omitted, NREP is 
set to the value of NPSU. (if both NPSU and NREP are 
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omitted, you are in some trouble.) In the current 
version, a maximum of 48 replicates can be created 
appropriately. If NREP is largey than this, the half- 
samples begin repeating the same design so that no new 
information is added. 


@ NDEP = <number of dependent variables>, this parameter must 
be omitted unless REGR is specified. It gives the 
number of dependent varialbes for which regressions 
are to be performed. The last NDEP variables in the 
variable list are taken as dependent variables. If 
REGR is specified, but NDEP is omitted, NDEP is set 
equal to l. 


VAR Statement 


VAR variables; 


The VAR statement lists the variables to be analyzed. All of the vari- 
ables listed must be numaric. If omitted, all numeric variables in the input 


dataset will be analyzed. 


STRATID Statement 


STRATID variable(s); 


The variables in the STRATID must be sufficient to uniquely identify 
the individual sampling stratum. The file must be sorted by the STRATID 
variables, otherwise execution will terminate. There should be exactly two 
PSUs within each sampling stratum. If only 1 is found, a "dummy" PSU with 
no cases is created for purposes of balancing the half-samples. If more than 
two PSU are found, the excess PSUs are ignored (after a warning message is 


printed). 


PSUID Statement 


PSUID variables; 


The variables in the PSUID list must be sufficient to identify the indi- 
vidual PSUs within each stratum. The file must be sorted by the variables 


in the PSUID list within ftath stratum. As indicated above, there must be 
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exactly two PSUs within each sampling stratum. If the sample did not follow 
a strict BRR design, it may be necessary to construct "pseudo strata" or 
“pseudo PSUs" in order to obtain appropriate results. Consult your local 


sampling statistician. 


WEIGHT Statement 


WEIGHT variable; 


The WEIGHT statement identifies the variable that is to be used as a 
weight for each case. One and only one weight variable must be specified. 
If the sample is unweighted, then a dummy weight variable must first be 
created. (If the sample is unweighted, do you really need a complicated 


algorithm like this?) 


DETAILS 


Missing Values 


BRRVAR's treatment of missing values depeads upon whether any multi- 
variate statistics are requested. If COV, CORR, or REGR are specified, 
listwise deletion of cases with missing data is always performed. If only 
univariate statistics are requested, BRRVAR deletes missing values sepa- 
rately for each variable, unless NOMISS is specified in the PROC BRRVAR 
statement. If NOMISS is specified, then listwise deletion of cases with 


any missing values is always performed. 


Missing values are not allowed for the weight variable. If any missing 
values are encountered, the program will terminate on the presumption that 


the weight variable was incorrectly specified. 


Printed Output 


The printed output from BRRVAR includes the following: 


@ a Summary Report indicating the number of cases processed 
and the number with missing values, 


e a design matrix (if requested) indicating which PSUs were 
included in each of the replicate half-samples, 


® all requested univariate statistics, including the overall 
estimate, the estimated sampling error of the overall esti- 
mate, and tne lower and upper 95% confidence bounds for the 
overall estimate, 


e if requested, an estimated correlation coeficient is printed 
for each pair of. variables in the variable list along with 
the estimated sampling error of each correlation coefficient 
and the lower and upper confidence bounds of each corre- 
lation coefficient, 


@ if requested, BRRVAR prints a relatively unformatted dump 
of the estimates generated from each of the replicate half- 
samples. The columns of this dump correspond to the half- 
samples and the rows to the statistics being estimated. The 
first column gives the overall estimate. 


EXAMPLE 


The attached example in APPENDIX A shows the JCL necessary to access 
BRRVAR at COMNET. In this example, a single stratum variable and a single 
PSU identifier are read. The input dataset contains five records for each 
PSU with three variables for each records The data are already sorted by 
PSU within stratum. 


PROC BRRVAR is called twice. In the first instance, simple univarite 
statistics are requested. The second request asks for both correlations and 
a regression run. In addition, the DESIGN and DUMP options are invoked to 
provide additional information. 


The printout resulting from the sample dataset is shown following the 
listing of the input records. 
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APPENDIX A 


Examples of an Input Dataset, SAS Statements and 
SAS Output 


EXAMPLE: Iuput Dataset 


00010 //X JOB (ED,AIR), ‘SAGE-WYOUNG’ ,PRTY=8, 

00020 // MSGCLASS=P ,NOTIFY=XEJQNY ,MSGLEVEL=(1,1) 
00030 /#ROUTE PRINT R207 

00040 // EXEC XSAS,OPTIONS=’S=727’ ,PRINT=P 

00050 //WORK . DD UNIT=SYSDA,SPACE=(CYL,(5,5)) 

00060 //SASLIB DD DSN=NCES.XEJQNY.SAGELIB, DISP=SHR 

00070 //SYSIN DD # 

0©O080 DATA TEMP: 

00090 INPUT STR 1 PSU 2 

QOQ100 4 3-4 Y S$-6 Z 7-8: 

00110 

00120 WT=4.0: 

00130 CARDS: 

00140 1A 
00150 1A 
00160 1A 
00170 1A 
00180 1A 
00190 1B 
00200 1B 
00210 1B 
00220 1B 
00230 1B 
00240 2A 
Q0250 2A 
00260 2A 
00270 2A 
00280 2A 
00290 2B 
00300 2B 
00310 2B 
00320 2B 
00330 2B 
00340 3A 
OO3S5S0 3A 
00360 3A 
00370 3A 
00380 3A 
aO390 3B 
00400 3B 
00410 3B 
00420 3B 
00430 3B 
00440 4A 
00450 4A 212 
00460 4A 314 
00470 4A 418 
00480 4A 516 
09490 4B 313 
00500 4B 412 
00510 4B 519 
00520 4B 621 
00530 4B 716 
00540 ;: 


WHF HSPWN PK DVOVONUNAHUOCKAPONUNUNAHWOKK DAWN Us 


KBNIOUPHUOUAUNKYNUOUPHUOUAWUNKNOUPHUOUAHONE 
UQPURNWDWONOUMYENWDPUNOUVUAHWDUSTNDWUVFPUNKWPUCNVAUN 
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00550 
00560 
003570 
00580 
00590 
00600 
00610 
00620 


EXAMPLE: SAS Statements 


PROC BRRYAR NREP=8 SUM MEAN STD; 


e 
? 


PROC 


GN DUMP: 


00630 
00640 
00650 
00660 
00670 
09680 


PROC 


VAR X ¥ Z: 
STRATID STR; 
PSUID PSU; 
WEIGHT WT: 


BRRVAR NPSU=9 NREP=8 MEAN SUM STD REGR NDEP=i DESI 


WAR X ¥ Zs 
STRATID STR; 
PSUID PSU; 
WEIGHT WT; 
CORR: 

WAR XY Z3 
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EXAMPLE: SAS Output 


STATISTICAL ANALYSIS SYSTEM 3:90 THURSDAY, MARCH 10, 1983 
PROC BRRVAR ~ BRR SAMPLING ERROR ESTIMATES 
SUMMARY INFORMATION 


40 OBSERVATIONS READ 
160 WEIGHTED OBSERVATIONS READ 


MISSING VALUES DELETED SEPARATELY FOR EACH VAR 


9 REPLICATES IN & PSU PAIR DESIGN . 


STATISTICAL ANALYSIS SYSTEM 3:30 THURSDAY, MARCH 10, 1983 2 
PROC BRRVAR - BRR SAMPLING ERROR ESTIMATES 
UNIVARIATE STATISTICS 


95% CONFIDENCE BOUNDS 


STATISTIC VARIABLE TD. N ESTIMATE STD. ERROR LOWER UPPER LABEL 
SUM x 160.0 640.000 80.000 +»=S «483.200 ©7298 ..800 
SUM y 160.0 1052.00 62.354 929,786 174.214 
SUM z 160.0 760.000 += 105.830 ©=»«-552,579 «947,427 
MEAN x 140.0 4.000 0.500 3,020 4,980 
MEAN y 160.0 6.575 0.390 5.811 7.939 
MEAN Z 160.0 4,750 0.661 9,454 6.046, 
STD xX 160.0 1.732 0.124 1.489 1.975 
STD Y 160.0 5.449 0.505 4,460 6.439 
STD 2 160.0 2,107 0.271 1,575 2.638 
16 
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EXAMPLE: SAS Output (cont.) 


STATISTICAL ANALYSIS SYSTEX 
PROC BRRVAR - BRR SAMPLING ERROR ESTIMATES 


HALF SANPLE DESIGN MATRIX 
Pol MATRIX 


1A 01010011 
18 LO101109 
ZA OO100E11 
28 11011000 
JA 01001110 
J8 LoLiooddt! 
4A OO01T1OL 
48 11100010 


7 PSUS SPECIFIED 
8 PGUS FOUND 
STATISTICAL ANALYSIS SYSTEX 
PROC BRRVGR - BRR SAMPLING ERROR ESTIMATES 
SUMMARY INFORMATION 
40 OBSERVATIONS READ 


160 HEIGHTED OBSERVATIONS READ 
Q REIGHTED OBS DELETED FOR MISSING VALUES 


9 REPLICATES IN 8 PSU PAIR DESIGN 
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4:30 THURSDAY, HARCH 10, 1983 


4:30 THURSDAY, HARCH 10, 1983 


d 


4 


EXAMPLE: SAS Output (cont.) 


STATIGOTICAL ANALYSIS SYSTEX 3:30 THURSDAY, HARCH 10, 1903 5 
PROC BRRVAR - BRR SAHPLING ERROR ESTIKATES 
UNIVARIATE STATISTICS 
954 CONFIDENCE BOLNDS 
STATISTIC VARIABLE © ATD. NsESTIMAVE STD, ERROR LOWER UPFER LABEL 


SUK X 160.0 = 640,000 80,00 483,200 794,800 
alld 100.0 1052.000 O2.35% 929,786 «174,214 
SUH t 180.0 780,000 105,830 $52,579 947,427 
MEAN X 160.0 4,000 1500 3,020 4,980 
NEAK Y 140.0 6.575 0.390 DAL 7,439 
MEAN 2 180.0 4,750 0,66 d,454 4,046 
STD X 140.0 1,732 0.124 1,489 1,975 
at) Y 140.0 JM? 0.505 4.400 6.499, 
aT 160.0 =. 2,107 0.271 1,97) 2,638 


STATISTICAL ANALYSIS SYSTEX 3:90 THURSDAY, MARCH 10, 1983 6 
PROC BRRVAR - BRR SAMPLING ERROR EST INATES 


CORRELATION ESTIMATES 
FIRST ©» SECOND = CORRELATION 95% CONFIDENCE BOUNDS 
VARIABLE VARIABLE = HTD, N ESTIMATE © STD. ERROR LOWER UPPER = PROB C0 
X 140.0 0,026 0,260 ~0,482 0.935 . 
i K 140.0 0.069 0.351 -0,618 0.756 
f Y 140.0 0,341 ),273 “0.194 0,877 
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EXAMPLE: SAS Output (cont.) 


STATIGVICAL ANALYSIS SYSTEH 4:30 THURSDAY, HARCH 10, 1983 
PROW BRRVAR - BRA SAMPLING ERROR ESTIMATES 
REGRESSION COEFFICIENTS 


FIRST =©=—»SECOND «REGRESSION 994 CONFIDENCE BOUNDS 
VARTABLE VARIABLE HTD, ESTIMATE © STD, ERROR —LOKER UPPER PROB £0 
2 160.0 0,039 0,410 “0.770 0,836 


X 
2 Y 180.0 0.942 0,284 “0.216 0,099 
MULT RSG 2 140.0 0,127 0,252 “0.366 0,621 


STATISTICAL ANALYSIS SYSTEH 3:0 THURSDAY, HARCH 10, 1993 
PROC BRRVAR ~ BRR SAMPLING ERROR ESTIMATES 
REPLICATE STATISTICS 
FIRST © SECOND §=—-REPLICATE § 99% CONFIDENCE BOUNDS 
VARIABLE VARIABLE WID. N ESTIMATE © STD, ERROR LONER UPPER PROB C0 


0.039 0.662 0,042 0.971 0.438 0,432 0,524 0,049 0,395 
0.342 0.048 -0,051 0.914 0,639 0,471 0.3! 0,260 0,800 
0.127 0.442 0,010 0.248 0.277 0.147 0.496 04179 0,590 
160,000 160.000 160.000 140.000 160.000 140.000 180.000 140,000 140,000 
L732 Ald 1.702 1.658 1,732 1,792 1,058 1.658 1,659 
Ded} 90783 5,886 5,925 4,928 4,766 4,924 4,030 5,052 
2107 2.000 2,000 1,792 2.449 1.732 2.000 1.732 2,294 
140,000 160:000 140,000 160.000 160,000 140.000 140.600 140.000 140,000 
0,026 -0.043 0.284 0,159 -0.281 0.915 -0.198 0.455 -0,145 
0,069 0.707 0,289 -0.£74 0,000 0,000 0.302 0,522 0,135 
0,341 0,048 0,047 0,307 0,588 06424 0.619 0.206 0,784 
640.000 800.000 440.000 720,000 640.000 440,000 54.000 540.000 540,000 
1092,0001112.0001128,0001072.000 974.0001072,0001032,0001088.000 934,000 
700.000 800.000 800.000 40.000 800.000 940.000 800.000 640.000 440,000 
0.000 0.000 0,000 0,000 0,000 0,000 0.000 0,000 0,000 
0,000 0.000 0,000 0,000 0.000 0.000 - 0.000 0.000 0,000 rrr 4a 
0,000 0.000 0,000 0,000 0,000 0,000 0.000 0,000 0,000 _ 
0,000 0.000 0,000 0.000 0.000 0.000 0.000 0.000 0.000 
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EXAMPLE: SAS Output (cont.) 


STATISTICAL ANALYSIS SYSTEM 3:30 THURSDAY, HARCH 10, 1999  § 


VAR TABLE N NEAN STD DEV SUH NINTHUM HAX TMU 
X 4) 4,00000000 175411004 140,00000000 1.00000000 7,00000000 
y 4) b.97500000 3691967459 263.00000000 1.000000 2100000000 
d 40 4,75000000 2419337340 190.00000000 1.00000000 9,00000000 


CORRELATION COEFFICIENTS / PROB > “R* UNDER HOsRHO=0 / N= 40 
X Y i 


X 1.00000 0.02649 0.06852 
0.0000 0.8711 0.4744 
Y 0.02649 1.00000 0.34138 
0.8711 . 0.0000 0,031 
1 ),06052 0.34138 1.00000 
0.6744 0.0911 0.0000 
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Program Listing 


nOdLO 
oOO0020 
oo030 
mood4ao 
OGOSS 
Naas 
oon7y 
OOOO 
Qog9 QO 
00160 
00210 
oOo1La0 
O0130 
OQ L4O 
OO250 
OoO1450 
O0170 
city 186 
M0190 
fara) 20 5) 
M9210 
oads20 
Y0200 
00240 
OoOZ50 
QHOZSO 
QOo270 
ONZBo 
90290 
C0300 
O90310 
o0320 
00330 
90340 
00350 
00360 
00370 
90380 
00390 
oo 400 
O0410 
00420 
00430 


00440 — 


00450 
00460 
00470 
00480 
00490 
OOOO 
00510 
NOS2O 
00530 
00540 
oOSS9 
N9549 
20370 
OUS8O 
00590 
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fH 


ff 


/#ROUTE PRINT 


//WHYSTEP 1 


(ED.AIR) . “D-WYOQUNG ’ SPRTY=8, 
MSGCLASS=P .NOTIFY=KEIGNY .MSGLEVELH(1_1) 
207 

XPLIXCL.SYSOUT=P. 


-C* 


// PSRM.PLI=’NM,ND.AL,NEST.X.NC(S) MAR(2.72.1) .NSYN(C) NOP, 

// P&kM.LKED=/MAP .XREF.LIST’. 

// COND.LKED=(4,LT,PLI) 

//PLI.S7SIN 

MAINsPROC OPTIONS(MAIN) : 

1: Soteaientontenteatentententententetentetetententontentntententententantentenntantentetenteteteatetentententeteteatetetentententetententen oy 
/* x / 
/* BRRYAR : S&mMPLING VARIANCE ESTIMATES USING tf 
/* BALANCED REPEATED REPLICATIONS (BRR}: H#/ 
/* +/ 
/* PROC BRRWAR < OPTIONS PARMS is %y 
i% YAR < WARIABLE LIST *: #/ 
fit STRATID <STRATUM IDENTIFIER(S3 af 
/% PSuUID “ PSU IDENTIFIER(S) =: tf 
/* NOTE: FILE MUST BE SORTED BY ¥/ 
/* STRATUM AND PSU IDENTIFIERS: #/ 
/* WEIGHT < WEIGHT VARIABLE =: %/ 
/* NOTE: WEIGHT VARIABLE LIST ™UST %/ 
/* CONTAIN EXACTLY 1 VARIABLE af 
/%* ADJBY < ADJ FILE IDENTIFIERS +: x/ 
/* - NOT YET IMPLEMENTED - tf 
/* MUST BE A SUBSET OF TRE PSU x / 
[1% IDENTIFIERS: / 
/* COMPUTE «< COMPUTED VARIABLES =: %/ 
/% - NOT YET FULLY IMPLEMENTED - +f 
/* %/ 
7% OPTIONS: */ 
/* “/ 
/# DESIGN-PRINT DESIGN MATRIX FOR HALF SAMPLES +/ 
/* DUMP-PRINT STATISTICS FOR EACH HALF SAMPLE +*+/ 
/* SUM-ESTIMATE SUM OF EACH VARTISBLE #/ 
/* MEAN-ESTIMATE MEAN OF EACH VARIABLE %/ 
/* STD -ESTIMATE STD GF EACH VARIABLE a/ 
/* COW -ESTIMATE COW MATRIX FOR THE VARTASBLES#/ 
/# CORR-ESTIMATE COR MATRIX FOR THE VARISFEBLES#/ 
/* REGR-PERFORM REGRESSION %/ 
/% NOTE: (1) LAST VARIABLE IN VAR LIST af 
/* IS THE DEPENDENT VARIABLE #/ 
/* (2) IF CORR IS SPECIFIED, */ 
/% STANDARDIZED BETAS ARE / 
/* ESTIMATED. OTHERNISE RAW x/ 
/* REGRESSION COEFFICIENTS ARE +/ 
/* SHOWN: x / 
/* #/ 
/* PARAMETERS: t/ 
/% %/ 
/% MPSU-TOTAL NUMBER OF PSU’S IN THE SAMPLE +f 
/* #/ 
/* NRED-—(QPTIONAL) NUMBER OF REPLICATE HALF- #/ 
/% SAMPLES TQ BE USED 4 / 
/% NOTE: IF OMITTED NREP If xf 
/* MIN(100.,NPSU/2) #/ 
/* NDEP NUMBER CF DEPENDENT VARIABLES IF at / 
/* REGRESSION IS SPECIFIED */ 
/* +f 


NOTE: THE LAST NDEP WARTABLES IN THE 


29600 | /% WAR LIST ARE TAKEN AS DEPENDENT #,y 


Oa6ia 4% VARISELES. EACH IS REGRESSED +f 
Na620 {> AGSINST THRE PREDICTOR (NON- s/f 
INEIQO /*# DEPENDENT?) VARIABLES. sf 
O9640 ft DEFAULT=1. s/f 
OB620 /i */ 
20660 ¢ EXTERNAL ROUTINES: 3/ 
00679 /* #/ 
COSEO /* SETPARM . #/ 
OOaoa 3 (4 ALLOC ¥/ 
DO7 OO /# PROCESS rd 
DoO71O + PRINTIT */ 
no 720 /# af. 
NO730 /* ard 
QO740 So er rr ee t/ 
DoOvoSO /* / 
POF GO /*# DECLARE LINK@GE TO EXTERNAL SAS ROUTINES af 
OO77O /* %f 
oO7BO DCL SASPLO ENTRY, 

G0790 UNUSED ENTRYCFIXED BIN(31)), 

PaBaO BYPASS ENTRY(FIKED BIN (31)), 

Ongia SASLOG ENTRYC(CFIXED BIN(31),FIXED BIN(31)) : 

noe2o /* #/ 
90830 f# DECLARE EXTERNAL ROUTINES.. #/ 
noB4a /%# *s t/ 
NoBsa DCL SETPARM ENTRY((#.#) FLOAT BIN(5S3)), 

ON8S0 ALLOC ENTRY((#.*#) FLOAT BIN(S3)), 

ROB? PRINTIT ENTRY((#.%) FLOAT BIN(S3)), 

Qo8g8o PROCESS ENTRY((#,*) FLOAT BIN(5S3)) : 

CO Bea /* tf 
o0900 /* DECLARE EXTERNAL PROGRAM VARS... #/ - 
NO91O /# i #/ 
00920 DCL (MEANFLG. 

90930 STDFLG. 

RO98O COVFLG, 

GO9Sa CORRFLG, 

oo9E60 REGRFLG. 

OO970 SUMFLG, 

vOIBO SSQFLG. 

00990 NOMSFLG. 

oL000 DESFLG, 

OLo10 DUMFLG, 

81020 OQUTFLG) BIT(i) EXTERNAL, 

Q1L9Ga CINPTR, 

01040 PSUPTR. 

OL050 PSUPTR2, 

O1060 STRPTR, 

OLO70 STRPTRE, 

o1080 OQUTPTR) POINTER EXTERNAL. 

01090 (NREP, 

91100 NY, 

011190 NSTAT, 

011220 NCOMP, 

Ci130 NDEFP, 

01146 NPSUID. 

01150 NSTRID. 

©1160 NADJID. 

O1179 NOES, 

O12B0 NMISS.NTOT) FIXED BIN(31) EXTERNAL, 

o1190 (WGTORBS .WGTMISS.NUMPSU) FLOAT BIN(53) EMTESNAL,. 


<) 
ERIC 
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NLSoo 
o1L210 
OLZ2C 
W1220 
01240 
O125e 
Gi2Z250 
O1270 
O1260 
OL29Q0 
fy e.% 
ni Goo 
01310 
02320 
012330 
22240 
Q1250 
02360 
91370 
01380 
Q13a90 
M1400 
o14i0 
01420 
1430 
01440 
01450 
O1460 
01470 
02480 
01490 
o1500 
021310 
o1520 
22530 
01540 
QO1550 
/ 
91560 
01570 
M1580 
91590 
014600 
91610 
91420 
01630 
01640 
01650 
21660 
01670 
91480 
914690 
o1700 
01710 
O1720 
O1730 
O1740 
O17350 
O1760 
91770 


re) 
ERIC 


ESTABL1.,1) FLOAT BIN(SS3) CONTROLLED : 


DOL I FIXED BIN(31) STATIC: 
/% af 
i* MESSAGES... #; 
/* ays 
DCL MSGi CHAR(8O) INITC’ WERSICN 2A‘’3>. 
MSGiS CHAR(8O) INIT(C’ PREPARED FOR N.C.E.S.’), 
MSGiB CHAR(80) INIT(’ BY THE STATISTICAL pile 
mSG1IC CHARCEO) -INIT< / GROUP IN EDUCATION ‘S4&G5 
MSG2 CHAR(80) INIT(’ SMERICAN INSTITUTES FoR RZSEARC 
MSG3 CHAR(80) INIT(’ P. OQ. BOX 1113”), 
MSG4 CHAR(8®) INIT¢’ PALO ALTO, CA 94392’), 
MSGS CHAR(80) INIT(C’ (415) 493-3550"), 
MSGF1 FIXED BIN(31) BASED(ADDR(MSG1)). 
MSGFIA FIXED BIN(31) BASED(CADDR(MSG&IA)), 
MSGFiIB FIXED BIN(31) BASED‘ADDR{(MSG1IEB)), 
MSGFiC FIXED BIN(S1) BASED( ADDR (MSGIC)). 
MSGF2 FIXED BIN(31) BASED(ADDR(MSG2)), 
MSGF3 FIXED BIN(31) BASED(CADDR(MSG3)), 
MSG&4 FIXED BIN(31) BASED(ADDR(MSG4)), 
MSGFS FIXED BIN(31) BASED(AEDR(MEGS) ): 
/* e/ 
/*% CREATE SAS ENVIRONMENT... THEN ROCK & ROLL.. #/ 
/%* */ 
CALL “SASPLO $ 
ALLOCATE ESTAB : 
CALL SETPARMC(ESTAB): 
CALL ALLOC(ESTAB): 
CALL PROCESS‘(ESTAB): 
COLL PRINTITC(ESTAB)s 
/* t/ 
/* PRINT MESSAGES AND KALT.. +/ 
/* x/ 
CALL BYPASS(1); 
IF I=0O THEN DO: 
CALL UNUSED(-1) =: . /% PRINT CORE MSG.. 
CALL SASLOG(MSGF1.,80) : 
CALL SASLOG(MSGF1A,80O) ; 
CALL SASLOG(MSGF1B.890) ; 
CALL SASLOG(MSGF1C.890) :; 
CALL SASLOG(MSGF2,80) ; 
CALL SASLOG(MSGF3.80) =: 
- CALL SASLOG(MSGF4,80) ;: 
CALL SASLOG(MSGFS.89) : 
CALL SASLOG(MSGF1.,80) : 
END: 
END MAIN : 
/*#MEM2%/ 
{* ‘ 
//LKED.SYSLMOD DD DSN=NCES.XEJQNY.SAGELIB,DISP=O0LD 


//LKED.SYSLIB DD DSN=SYS2.PLIXLINK ,DISP=SHR 


// DD DSN=SYS2.PLIXBASE,DISP=SHR 

‘es DD DSN=SYS2.PLIXTASK.DIS®=SHR 

cf DD DSN=NCES.XEJGNY.SAGELIB, DISP=SHR 
/f DD DSN=SYS2.SAS.V796.SUBLIB.DISP=SHR 
//LKED.SYSIN DD +# 


ENTRY ENTRY 
INCLUDE S'YSLIB(SASPLO,PROCESS,ALLOC.SETPARM) 
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01780 NAME BRRYVARZ(R) 
QL790 /* 
END GF DATA 


9 


a1 S/H -OB (ED.AIR). /O-WYOCUNG’ LPRTY=8, 

woz f// MSGCLASS=P ,NOTIFY=SHXEJQNY MSGLEWSta(1 4) 
O30 /SRQUTE PRINT R2EO7 

10040 //WYYSTEPL EXEC XPLIXCL.SYSOUT=F, 

0050 // PARM.PLI=’NM_.ND.A.NEST.X.NC(S) MARC2. 72.1) .NSYN(C) NO” OFC’. 
MIS // PARM.LKED=’MAP.KREF.LIST.,NCAL’. 
1070 // COND.LKED=(4_,LT.PLTI) 

19080 //PLI.SYSIN DD + 

19090 ALLOC:PROCCESTAR) 


19100 ea ne a ait aoe eaten oleate t/ 

19210 (* . a/ 

19120 /# ALLOC-DETERMINES AMQUNT OF CORE PEUGIFPED FOR %/ 

1125 /* NASSVAR DATA STRUCTURES.. #/ 

101406 /* x/ 

10186 st ALLOCATES FPERUIRED CORE =/ 

26269 /+* */ 

i270 /% ZEROS CORES ALLOCATED e/ 

LEO i #/ 

10190 f* EXTERNAL FUNCTION ‘ALLOTRN’ [9 IMVoNeD x/ 

16200 /* TO DETERMINE IF COMPYAR AND OUTVAR af 

JO2Z720 /* STSTEMENTS ARE USED... “ALLCTRN ’ %/ 

0220 / TAKES CARE OF ALLOCATING CORE FOR THESE af 

JO230 /% STRUCTURES... ‘ALLCTRN’ RETURNS /’1°R IF x / 

O240—— /* THE COMPVAR/OUTYVAR CONSTRUCTION I uUSED,, #/ 

290250 /* #/ 

10260 Pm er rrr rrr t / 

W270 /* : #/ 

IN280 /# “ARE EXTERNAL LINKAGE TO SAS.. 4 / 

10290 /* 3 H/ 

30300 DCL NLJA ENTRYCFIXED BIN(31)) 

20310 RETURNS(FIXED BIN(31)). 

30320 GETMEM ENTRY(FIXED BIN(31).FIXED BIN(31)3.°TR. 

IO330 FIXED BIN(31)), 

ING4d ZERO ENTRY(PTR,FIXED SIN(31)), 

INGO MEMERR ENTRY(FIXKED BIN(31)). 

30360 LODOUT ENTRY: 

20370 4% tf 

IN3B0 /* EXTERNAL VARTABLES... %/ 

29390 /* f/f 

20400 DCL (INPTR,.PSUPTR,PSUPTR2.STRPTR.STRPTR2.OUTPTR? 

JO41O PTR EXTERNAL. 

20420 (COVFLG .CORRFELG.REGRFLG .OQUTFLG.SSCFLG.NOMSFLG) 

20430 BIT(1) EXTERNAL. 

30440 (MEANFLG,STDFLG.SUMFLG) BIT(1) EXTERNAL. 

399450 - (NSTAT.NREP.NYV1,NPSUID,.NADJID.NSTRID) 

20460 FIXED BIN(31) EXTERNAL, 

29470 ESTAR(#,#) FLOAT BIN(S2) CONTROLLED : 

29480 /* : ‘ / 

00490 /% LOCALS.. #/ 

20500 /* #/ 

9O510 DCL (NY NUE. NREPL.NEED.GOT) FIXED BIN(31), 

IOS20 WORKYVEC(1) FLOAT BIN(S3) BASED(INPTR), 

00530 PTR PTR. 

290540 NULL BUILTIN: 

JO550 /* %/ 

00560 7* FREE ESTIMATE TAB... PATCH TQ RUN UNDER THE : %/ 

J03570 /*# OPTIMIZER WHICH REQUIRES THAT A CONTROLLED ef 

00580 /# STRUCTURE MUST BE ALLOCATED IN ORDER TO BE + / 
g 18870 /* PASSED... (NOT THE “OF WITH THE ‘F’ COMPILER... x/ 

ont 


9600 
QOO610 
MIG2O 
99620 
00640 
OO6S0 
mN66% 
NQ670 
C°e689 
ON69O, 
ga700 
Qa710 
cad720 
C9720 
On740 
O07350 
09760 
Q0775 
00780 
C0790 
ANPOC 
ooela 
COB2ZO 
00830 
QOCE40 
Oossoa 
90860 
00870 
cogeo 
OOBFO 
QOF1LOG 
gO929 
09930 
00940 
09950 
0960 
00970 
QO°BO 
M0990 
91000 
01010 
01026 
01020 
91040 
Q1050 
01060 
01070 
01080 
01090 
O1100 
01110 
01120 
¢11230 
01140 
01150 
01160 
O1170 
01180 
01190 


<) 
ERIC 


/*® ST THIS POINT WE HAWE NO IDEA NHAT DIMENSIONS 

/* THIS TABLE MUST HAWE.... WILL ALLOCATE IT LATER 
f/%* $= IN THIS MODULE.... 

/% 

FREES ESTAB : 

/% LOAD SAS OUTPUT MODS If QUTPUT SPECIFIED BEFORE 
f* ALLOCATING ANY CORE(QUTMODS NEED ABOUT t#..) 

S% 

QUTFLE&=/0’B: © /# NOT YET IMPLEMENTED */ 

IF OUTFLG THEN CALL LODOUT : 


/* 

/* NEED=DOUBLE WORDS REGUIRED FOR TEMPORY 
/* STORAGE 

/# 

/* NWI-FOR INPUT VARS 

/* NPSUID-FOR PSU ID VYARS 

(% NPSUID-FOR PRIOR FSU ID 

/% NSTRID-FOR STRAT ID WARS 

/* NSTRID-FCR PRIOR STRAT ID 

/* 


NEED=B4ANV14+164#NPSUID+1G64#NADJID+16#NV1: 
CALL GETMEM(NEED,. GOT, INPTR.O): 
IF NEED < GOT THEN DO: 

CALL MEMERR(O): 

STOP: 
END: 
PSUPTR=ADDR(WORKVEC(NUL4+1))5 
PSUPTR2=ADDR(WORKVEC(NVY14+NPSUID+1))s 
STRPTR=ADDR(CWORKVEC(NVI+24NPSUID+1)): 
STRPTRZ=ADDRIWORKVEC (NYVI+24NPSUID+NSTRID+1) ):° 


/* 

/* ZERO REGION RETURNED... 

/* 

CALL ZERO(INPTR.NEED) : 

/* 

/* NSTAT=DOUBLE HORDS REQUIRED FOR EACH 

/* REPLICATE HALF SAMPLE: 

/* NULSENU+1 

/* NUZENVANVI/2 

/* . NY1 -— FOR RECORD COUNTS (CIF \NOMISS) 
/* NY¥V1 - FOR WEIGHTED SUMS (IF \NOMISS) 
/* N¥V1I+NCOMP - FOR SUMS/MEANS 

/* NW - FOR SSQS/STDS 

/* NU2 - FOR SSCP/COV/CORR 

/* 

NWENUI-1s 


IF COVFLG * CORRFLG * REGRFLG THEN DO: 
SSGFLG='17Hs 
NYZ=NVLIENYU/2 3 

END: 

ELSE DO: 

SSQFLG='0°Es: 
NW2=03 


END? 


NSTAT=24#NV14+NV2: 

IF NOMSFLG THEN NSTAT=NSTAT+4: 
ELSE NSTAT=NSTAT+24#NV15 

OUTPTR=NULL: 

/% OUTPUT NOT SUPPORTED 
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%/ 


tf 
+ / 
%/ 


*/ 
a] 
xs 
*/ 
%/ 
*/ 
*/ 
#/ 
*/ 
#/ 
/ 


tt / 


012900 
01210 
921220 
C1230 
giZ4aa 
O1250 
01260 
O1270 
01280 
01290 
01300 


/* xf 
/+ ALLOCATE 2 DIM TS&EBLE TO HOLD ESTIMATES... a; 
/* af 


NREPLSNRES+1: 
NEETSSH#NREPL#NSTAT: 
CALL GETMEM(-1,GOT,PTR,O) : 
IF NEED > GOT THEN DO : 
CALL MEMERR(NEESD-GOT) =: 
STQP s 
END : 
ALLOCATE ESTAR(NREP1.NSTAT): 


Qi3190 END ALLOC : 


Q1320 /*#MEM4#/ 


O1IIE /# 


02340 //LKED.SYSLMOD DD DSN=NCES.XEJQNY.SAGEL IE DISP=SHR 
01350 //LKED.SYSLIB DD 


01360 // 
O2G70 // 
o1280 // 
02390 // 
o1da0 // 
O2r410 // 
O2420 // 
O1430 // 


DD 

DD DSN=NCES.XEJQNY.SAGELTE.DISP=SHR 

DD DSN=SYS2.CAS.V7946.SUBLIB.DISP=SHR 
DD DSN=SYS2.SAS.R795.GLTBRARY .DISP=SHR 
DD DSN=SYS2.PLI? SASE.DISP=SHR 


O1440 //LKED.SYSIN DD # 
O1450 NAME ALLOC(R) 


M1450 /%# 
END OF DATA 


re) 
ERIC 
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re) 
ERIC 


QOG1LO 
QOGo2z9 
eaa2ng 
90040 
GHG 
AgGso 
OHO7O 
Oooggo 
ONN90 
20100 
oo110 
OoL20 
OO2130 
oo140 
OyLeo 
00150 
QO170 
Oo180 
Q0190 
go20e 
on2ia 
Mg220 
WO230 
0240 
OS5O 
NO2EQ 
OO270 
QOO2BO0 
OG290 
gosaa 
o0o310 
00S20 
0O330 
O0340 
OOS 
00360 
00970 
OOSEO 
QN390 
Q0400 
a a 

oo410 
00420 
OO43a 
00446 
O0450 
OC460 
09470 
00480 
OG499 
QOSG0 
OoS1G 
00520 
OOSSO 
OO0S4a 
QOS50 
00560 
00570 
OOSBO 


{/% 
// 


//WVYYST=P1 EX 


// PARM.PLI=’NM,ND.,A,NEST.X.NC(S) MAR(Z, 72.1) .MS¥N(C) NOP LCE oe, 
// PARAM. LAED= ‘MAP KRESS. LIST. NCOAL’, 
if COND.LKED=(4,LT,PL2) 
//FLI.SYSIN DD # 
SETPARM:PROC(ESTAB) ;: 
| oe -- -- - - $e +--+ ef 
/* 4/ 
/*- SETPARM-DISCOVERS OROCEDURE OPTIONS AND PARMS #/ 
/%* VALIDATES SAMPLE DESIGN... a/ 
/* af 
[ fo men - $$$ $$$ 5 = 5 = - +--+ -- = === - */ 
/* x/ 
/%* DECLARE EXTERNAL LINKAGES... 4/ 
/* x/ 
DCL IOPT ENTRYCFIKED BIN(31)) 
RETURNS(FIXED BIN(31)), 
PARM ENTRYCFIXED BIN(G1)) RETURNS( FLOAT BIN(SS)? 
ERROR ENTRYC(FIXED BIN(31). *,#) FLOAT BIN(S3)) 
OBSERR ENTRY, 
NOVAR ENTRY(FIXED BIN(31)) RETURNS(FIKED BIN(3i3) 3: 
/* a/ 
/%* DECLARE EXTERNAL PROGRAM VARIABLES. uf 
/ * #/ 
DCL (NREP,NW1 .NDEP.NPSUID,NSTRID,NADJID,NCOMP) 
FIXED BIN(31) EXTERNAL, 
NUMPSU FLOAT BIN(53) EXTERNAL, 
ESTAB(#,#) FLOAT BIN(S3) CONTROLLED. 
(SSQFLG,CORRFLG.MEANFLG STDFLG,COVFLG,REGRFLG.SUMFLG, 
DESFLG,DUMFLG,NOMSFLG) BIT(1) EXTERNAL: 
/* af 
/%* ERROR MESSAGES... %/ 
/* x/ 
DCL ERR1I CHAR(809) INIT( : 
‘ ERROR: PSU IDENTIFIER(S) NOT SPECIFIED.’), 
ERR2 CHAR(80) INIT( 
‘ ERROR: ADJ BY LIST NOT COMPATIBLE WITH PSU ID LIST, 
ERRS CHAR(8O) INIT( 
’ ERROR: NO WEIGHT GARIABLES SPECIFIED. ’) 
ERR4 CHAR(89O) INIT( 
’ ERROR: NUMPSU PARM MISSING OR INVALID.’), 
ERRS CHAR(8O) INIT ( 
’ ERROR: MORE THAN 64 REPLICATES SPECIFIED.’ ). 
ERRS CHAR(80) INIT( 
’ ERROR: STRATUM IDENTIFIER(S) NOT SPECIFIED.’), 
ERR7 CHAR(80) INIT ‘ 
’ ERROR: INVALID SPECIFICATION FOR NDEP.’), 
ERRIF FIXED BIN(31) BASED(ADDR(ERRi)), 
ERR2F FIXED BIN(31) BASED(ADDR(ERRS)), 
ERRSF FIXED BIN(31) BASED(ADDR(ERRS)), 
ERR4F FIXED BIN(G1) BASED(ADDR(ERR4)). 
ERRSF FIXED BIN(31) BASED(ADDR(ERRS)), 
. ERR6F FIXED BIN(31) BASED(ADDR(ERRG)} , 
ERS7F FIXED BIN(31) BASED(ADDR(ERR7))s 
/%* #/ 


MSGCLASS=F ,NOTIF 
/#ROUTE PRINT R207 


ea 


JGB (ED,AIR), ‘D-NYOUNG’ JPRTY=8, 


YeXEsQNY ,MSGLEWVEL=(1.2) 


C XPLIXCL,SYSOUT=P., 


OoS90 
OUs6OC 
19610 
90620 
(ISSO 
Q0640 
M9650 
90640 
ler wae 
02680 
C0690 
gdg700 
co710 
oo720 
00730 
00740 
oo750 
Ua760 
OO770 
00780 
coa79Q 
oGg00 
aog10 
Oca zo 
GOB3a 
90840 
COBESO 
ocgéa 
00870 
o7B880 
90290 
09900 
00910 
O0920 
90930 
00940 
00950 
90960 
00970 
q0o9°80 
00990 
C2000 
01010 
01020 
01030 
O1040 
01050 
01060 
01070 
O1080 
91090 
01100 
01110 
01120 
QL1LSO 
212140 


OLI50 - 


31160 
01170 
01180 


re) 
ERIC 


/% DISCOVER PROCEDURE OPTIONS.. xf 
/% af 
IF I0PT(1)=+1 THEN MEANFLG='1°B 

ELSE MEANFLG=/0°R ;: 

IF YOPT(2)=1 THEN STDFLG='1°8 

ELSE STDFLG='0°’E : 

TF TOPT(S)=1 THEN COVFLG=/1°s 

ELSE COVFLG='0°R ; 

IF IOPT(4)=1 THEN CORRFLG=‘/1°8 : 

ELSE CORRFLG=/0°'R ; 

IF IOPT(S)=1 THEN REGRFLG=/1°' 

ELSE REGRFLG='0°'R : 

IF IOPT(6)=1 THEN SUMFPLG=/1°E: 

ELSE SUMFLG='0O’R; 

IF IOPT(7)=1 THEN DESFLG=/1'8: 

ELSE DESFLG='0'R: 

IF IOPT(S)=1 THEN DUMFLG='1°8 : 
ELSE DUMFLG='0°R : 

IF IGP7(9)=1 THEN NOMSFLG=/1°H: 

ELSE NOMSFLG='0°H; 


/* */ 
/* INSURE COMPATIBILITY OF OPTIONS & SET DEFAULTS ¥/ 
/* tf 


IF \(MEANFLG*STDFLG*COVFLG“CORRFLG*REGRFLG) 
THEN SUMFLG=‘1/‘H: 
Ir COVFLG°CORRFLG“REGRFLG THEN SSQFLG='1°‘H: 
ELSE SSGQFLG='0O'8: 
IF SSGFLG THEN NOMSFLG='1'H: 
/* DETERMINE NUMEER OF VARS ON EACH LIST #/ 
NPSUID=NOVAR(3); 
IF NPSUID < 1 THEN DO: 


CALL ERROR(ERRIF,.ESTAE): 
END: 
NADJID=NOVAR(4): 
IF NADJID > NPSUID * NADJID > O 


THEN DO: /* RATIO ADJUSTMENTS #/ 
/* NOT YET HERE + / 
CALL ERROR(ERR2F.ESTAB): 
END: 


NSTRID=NOWAR(S): 
IF NSTRID < 1 THEN DO: 
CALL ERROR(ERR6F,ESTAB): 


END: 
/* a/ 
/%* DETERMINE NUMBER WARS ON WEIGHT STATEMENT.. + / 
/* x/ 


IF NOVAR(6) \= 1 THEN DO ; 
CALL ERROR(ERRSF,ESTAB) ;: 
END 


/* DETERMINE NO OF COMPUTED VARS t/ 
NCOMP=0: 

NCOMP=NOWUAR(C7): 

/% #/ 
/* SET NUMPSU %/ 
/* %/ 


NUMPSU=0L00: 
NUMPSU=PARM(1): 
/* #/ 
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01190 
Q1200 
QO1210 
O1220 
O1230 
01240 
o1250 
O1260 
G1270 
o12ec 
01290 
Q1200 
019310 
O2020 
* 91930 
Q1040 
OLOS0 
M12é6c 
91379 
O1380 
O1390 
91400 
01410 
01420 
G1430 
O146a 
O1450 
O1440 
01470 
O1460 


/% CHECK NREP 

/* 

NREP=0.0: 

NREP=PARM: 2): 

IF NREP < 41 THEN NREP=(NUMPSU/2+1): 
ELSE NREP=NREP +1: 

NULENOVYAR( LI +15 

/% CHECK NDEP 

NDEP=0: 

NDEP=PARM( 3): 


IF REGRFtLG THEN IF NDEP « 1 THEN NDEP=1: 
IF NDEP >= NQWAR( 1) * (\REGRFLG & NDEP } O) 


CALL ERROR(ERR7F .ESTAB): 

END: , 
END SETPARM : 
/AMEMNQH/ 

/* 


/f/LKED.SYSLIE DD 


//LKED.SYSIN DD # 
NAME SETPARMC(R) 
/# 


END OF DATA 


re) 
1 oeeeneneeenee 
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ere nner td erent nrtnrn epn 


//LKED.SYSLMOD DD DSN=NCES.XEJGNY.SAGELIE, DISP=0LD 


At DD 

// DD 

// DD 

// D2 

// DD DSN=NCES.XEJGNY.SAGELIB, DISP=SHR 

// DD DSN=SYS2.SAS.V796.SUBLIE, DISP=SHR 
f/ DD DSN=SYS2.SAS.R79S.GLIBRARY.OfSP=SHR 
/ DD DSN=SYS2.PLIXBASE,DISP=SHR 


THEN 


D 


QO: 


x/ 
%/ 


NQa1o 
WoOLO 
00030 
00040 
oo0saC 
o3060 
OOU7O 
so0aBO 
WIet@h- 28) 
OOLd0 
QOL 
oc120 
NO1S0 
oo14d0 
OO1S0 
CO160 
oO170 
001980 
oOoLg9a 
09200 
OOZi0 
QO220 
COo2Zoo 
O0240 
O02Z80 
09250 
00270 
o02ZB0 
BOZFO 
GOB 
SOSL0 
e0320 
00330 
00340 
OO350 


_ 00360 


re) 
ERIC 


00370 
00380 
00390 
90400 
00410 
09420 
00420 
00440 
N0450 
90460 
00470 
90480 
00490 
90500 
Nes1o 
00320 
00530 
Nag540 
OOSEO 
00560 
oos7a 
00380 
OOS90 


ars JOB 


ff 


(ED.,AIR). ’D-WYOUNG’ ,PRTY=8, 
MSGCLASS=P ,NOTIFY=xXEVJQNY .MSGLEWEL=(1,1) 


/*#ROUTE PRINT R207 

//WYYSTEP? EXEC 4PLIXCL.SYSOUT=P. 
// PARM.PLIS=’NM.ND,A,NEST,X,NC(S) MAR(2,72.1) .NSYN(C) NOP, OFC’. 
// PARM.LKED=‘/MGP.xXREF,LIST,NCAL’, 
// COND.LKED=(4,LT.PLI) 


f/PLI.SYSIN 


DD + 


PROCESS: PROC(ESTARB) : 


PO rer er ne et f 
/% tf 
/* PROCESS-COMPUTES ESTIMATES OF CHARACTERISTICS a/ 
/%* SPECIFIED FOR EACH REPLICATE IN THE wf 
/* SAMPLE DeESIGN.. x 
/* es 
/* WRITES REPLICATE LEVEL RECORDS TO OUTPUT #7 
/% DS IF THE OUTPUT OPTIGN [S$ IN EFFECT uf 
/* #/ 
/* COMPILES AND EXECUTES PROESAMmMING 2%, 
/*% STATEMENTS FOLLOWING THE PROC IF */ 
(%* THE OQUTVAR/COMPYAR CONSTRUCTION a} 
/* IS WSED.. ; #/ 
/* af 
/* NATIONAL ESTIMATES ARE STORED IN #/ 
/* ARRAY ESTAB(NREPS .VARNO).. ef 
J af 
[ Romeo -- --  - - 5-5 5 #/ 
/* _&/ 
/* EXTERNALS.. tf 
/* ue 
DCL (OUTFLG, QUTYFLG,SSGFLG.NOMSFLG,DESFLG.DUMFLG) 

BIT(1) EXTERNAL. 

(OUTPTR. INPTR,STRPTRISTRPTR2, 

PSUPTR.PSUPTRZ) PTR EXTERNAL, 

(NOBS.NW1 JNSTAT,NREP .NPSUTD.NMISS.NSTRID) 

FIXED BIN(31) EXTERNAL, 

(WGTOBS, NUMPSU.WGTMISS) FLOAT BIN(S2) EXTERNAL 

/* #/ 
/* SAS LINKAGE... u/ 
/# #/ 
DCL INPUT ENTRY RETURNS(FIXED BIN(21)), 

WARX ENTRY(FIXED BIN(31).FLOAT BIN(S2)), 

OBSPTR ENTRY RETURNS(PTR), 

SETDSN ENTRY(FIXED BIN(31)?). 

PUTOUT ENTRY(FIXED BIN(3!),PTR). 

FRUF ENTRY, 

BYPASS ENTRY(FIXED BIN(31)), 

CLSOUT ENTRY. 

IOPT ENTRY(FIXED BIN(31)) 

RETURNS(FIXED BIN(31)?, 

LINSIZ ENTRY(FIXED BIN(31)). 

STITLE ENTRY(FIXED BIN(31),FIXED BIN(31)). 

NAMEY ENTRY(FIXED BIN(31),FIKED BIN(31), 

FIXED BIN(15).). 

ERROR ENTRY(FIXED BIN(31).(#.#) FLOAT BIN(S3)), 

MISSING ENTRY(FLOAT BIN(S3)) RETURNS(BIT(1)): 
/* 4/ 
/* LOCALS.. %/ 
/* s/f 
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re) 
ERIC 


WOEOO 
Cagle 
coe620 
mO63o 
oné4a 
Qd6a0 
99660 
00579 
oogé6so 
00690 
00700 
Oo710 
00720 
00730 
voO740 
dO750 
COTEC 
00770 
@TIC: 
OO7RCa 
OO790 
oosco 
90810 
OOBZO 
an 
CORSO 
Oo’. 
Owe4o 
YOBIO 
48). 
O02 60 
11). 


008709. 


yO8BO 
ogeg90 
20900 
OOF19 
00920 
00930 
00940 
OO9S0 
00960 
00970 
00980 
00990 
01000 
o1010 
o1020 
01030 
01040 
01050 
o1060 
Q10709 
01080 
91099 
01100 
O1140 
01120 
02130 
01140 


DCL 


DCL 


DCL 


/*® 
/% 
/* 
DCL 


DCL 


ESTARB(#.#) 
WGT FLOAT 


FLOAT BIN(S3) CONTROLLED, 
BIN(S3) STATIC, 


NPSU FLOAT BIN(SS3) STATIC INIT(O). 


INHOLD(1> 
EYPSU(1) 
BYPSUZ(1) 
BYCHPSUC1) 
BYCHPS2(1) 


FLOAT BIN(S3) BASEDCINPTR). 

FLOAT BIN(S3) BASED(PSUPTR). 

FLOAT BIN(S2) BASED(PSUPTR2>, 
CHARS) BASEDCPSUPTR). 
CHAR (8) BASED(PSUPTR2), 


BYSTR(1) FLOAT BIN (53) BASED(STRPTR), 


BYSTR2(1) 


FLOAT BIN(S3) BASED(STRPTR2) , 


BYCHAR2(64) BIT(1) STATIC, 


BYCHSTR(1) 
BYCHST2(1).) 
BYCHAR( 64) 


CHAR(G) BASED(STRPTR), 
CHAR(8) BASED(STRPTR2), 
BIT(1) STATIC, 


OUTYVEC FIKED BIN(31) BASED(OQUTPTR), 


(COPTR,IPTR) 


PTR, 


(ADDR, SUBSTR) BUILTIN, 


CIR, IV. K NY KY, IRI, JU.L,LNSZ. CENTER) FIXED 


DESIGN(16) 


CHAR(8) STATIC INIT<( 


/%® B%/ %11101900K", 


{*#16%/ 74 


‘40 
/824R/ 4111121010", (1190901100. 419010000K’, 
/#O28%/ (MOGO1TO1O’. fA 


/#48%/ 4 


2d 


MXREPS(S) 


REPIND(S) 


1110101°, OLdOaK ”, 


101100, ‘OL111100%, “Limi9dd 


SAT DOL!. “SPOC TOTO", *21L0C1O0*. *2102109 


O10L100/. “OOLOOOaK! 3s 
FIXED BIN(31) STATIC INIT(S. 16, 24, 32, 


FINED BIN(31) STATIC INIT(1. .2. 4, 7. 


(X, WT, XWT.REPLID) FLOAT BIN(S3) STATIC. 


DESVEC 
(STRATNO,. 
{LASTPSU, 


CHAR(64) BASED(DESPTR). 
ISINDX, IRINDX} FIMED BIN(G1) STATIC, 
TESTPSU) CHAR(1) STATIC, 


DESERR BIT(1) STATIC, 


DESPTR 
FLAT_EST(1 
NAMESTR, 
NTYPE 
NPOS 
NLNG 
NVARO 
NNAME 
NLABEL 
NFORM 
NIFORM 
NFL 

NFD 

NF 

NJUST 


Hthn Nn Nh Nh hoe 


ERROR MES 


ERR1 CHAR ¢ 
‘ ERROR: 
ERRIF FIXKE 
ERR2 CHAR 
‘ ERROR: 


POINTER, 
» FLOAT BIN(S3), 


FIXED BIN(1S5), 

FIXED BIN(15), 

FIXED BIN(15), 

FIXED BIN(15), 

CHAR(&), 

CHAR( 40), 

CHAR(S), 

CHAR(8), 

FIXED BIN(15), 

FIXED BIN(15), 

FIXED BIN(15), 

FIXED BIN(15) 

%/ 
SAGES.. te/ 
*/ 

80) INIT ( 

MISSING VALUE FOR WEIGHT INVALID. /?}. 
D BIN( 21) BASED(ADDR(ERRI)): 

(80) STATIC INIT¢ 
FILE NOT SORTED BY STRAT & PSU ID WARIAELE 
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Su’). 
01250 
91160 
o1170 
01280 
Y 
01190 
O1200 
OL1210 
02220 
O125% 
GL240 
o1250 
02260 
OLe7o 
Oi 286 
GL290 
01300 
91310 
O2920 
O1290 
o13940 
O1L350 
91960 
01370 
01380 
01390 
OL400 
o1410 
01420 
01430 
01440 
01450 
01440 
01470 
01480 
01490 
C1500 
01510 
01520 
01530 
C1540 
O1550 
01560 
O1570 
M1580 
OLSE0 
01600 
01610 
01620 
01630 
21640 
01650 
01660 
1670 
0168¢ 
OL690 
01700 
01710 
01720 


re) 
PRiGs as 


SETDES: 


ERRZF 


FIXED BIN(G1) BASED(ADDR(ERR2) 


CN ENDPAGE(SYSPRINT) BEGIN 3: 


CALL ST 


IPLetOvL> 


at 


PUT EDIT( ‘PROC SRRYAR ~- BRR SAMPLING ERROR ESTIMATES’ 


PUT SKI 
END; 


/* SET LI 
/* 

CALL LINS 
CENTER=(LN 
IF IQPT(SO 


/* IF QU 
/% 
IF OUTFLG 
/% 
/* IN 
/* 
DO Iv 
CAL 
IF 


ELS 
END: 
pO ive 
CALL 
IF N 


END: 
LASTPSU 
NORS=0: 
DO IR=1 

po I 

E 

END 

END : 

NVENU 1 = 
IR1=NRE 
STRATNO=0; 
LASTPSU=/1 


DO IR = i 
IF NREP 

ENDs 

IR=S: 


IVY=REP IND ¢ 
DESPTR=ADD 
ISINDX=1: 
IF DESFLG 


(SKIP(1) .COL(CENTER).A) : 
P (2) 3 


MNESIZE AND CENTER 


TZ(LNSZ): 


SZ-35)/2: 
) = 0 THEN CENTER=1 


TPUT SPECIFIED DEFINE QUTFUT DS 


THEN CALL PREFOUT 
ITIALIZE VALUES 


1 TO NPSUiID: 

L NAMEW(S, IV NTYPE ,NAMESTHR) : 
NTYPE=2 TREN 
BYCHARCIY)=/1°R: 

EF BYCHARCIV)=/0°'R: 


TO NSTRID: 
NAMEV(S,IYV.NTYPE,NAMESTR); 
TYPE=2 THEN BYCHARZ( IV} =/1°E: 

ELSE BYCHAR2Z(IV)='9O'B: 


=/1’Es 

WGTORS=0;3 

TO NREP+1 $ 

Y=l TO NSTAT: 
STARCIR, 1¥)=0.0 + 


« 
Y 


1; 

P+1: 

TO 3S: 

<= MXREPS(IRI+2 THEN GO TO SE 


IR): 
R(DESIGN(IW)): 


THEH DO: 


SIGNAL ENDPAGE(SYSPRINT): 
PUT EDIT( ’HALF SAMPLE DESIGN MATRIX’, 


END; 


PUT SKIP; 


/# 


(COL(CENTER) ,A.SKIP (2) .A.CGL¢ 
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“PSL 
20) 


a / 


é 


& YARS.. / 


~/ MATRIX’) 


Ad 


*/ 


4 


m) 
ERIC 


01730 
01740 
OL7E0 
QO176o 
01770 
OL7sa 
OL7F0 
o2rB800 
01810 
01829 
OLB8ga 
O1849 
OiBsa 
91850 
01970 
o18g0 
01890 
M1900 
Q1910 
01920 
91930 
QL94AS 
01950 
W1960 
01970 
C1986 
01990 
Us: 

O2000 
STAB) 
02010 
02029 


02030 
TAB): 
92040 
02050 
02060 
92070 
QZ2080 
O2090 
02100 
02110 
02120 
02136 
02140 
02150 
02160 
02170 
G2180 
92190 
02200 
02210 
o2220 
02230 
02240 
02250 
02260 
02279 
022680 


fs STEP THRU OBS & COMPUTE ESTIMATES AT THE e/ 
/* REPLICATE LEVEL... t/ 
ft uf 


DO WHILE( INPUT=0) + 
CALL VARM(S. BYSTRIL)): 
CALL VARK(3. BYPSU(i)): 


IF NCBS = 0 THEN GO TQ NEKSTR: 
f* CHECK FOR NEXT STRATUM * ef 
DO f¥=1 TO NSTRID: 
IF BYCHAR2( IY) THEN DQ: 
IF BYCHSTR(IV) » BYCHST2( IW) 
THEN GO TO NENSTR: 
ELSE IF BYCHSTR(IW) « BYCHST2( IW) 
THEN CALL ERROR(CERR2ZF ,ESTAB): 


ELSE DO: 
IF BYSTR(IV) > BYSTR2( IW) 
THEN GO TO NEWSTR: 
ELSE IF BYSTR(IV) < BYSTRACTW) 
THEN CALL ERROR(EPRES ESTAR): 


{= CHECK FCR NEXT PSU #/ 
DO IY = i TO NPSUID;: 
SF BYCHARCIV) THEN DO: 
IF BYCHPSUC IV) > BYCHPS2(IW) THEN GO To NEWPS 


ELSE IF BYCHPSUCIVI<BYCHPS2( IV) THEN CALL ERROR(ERR2F LE 


GO TO 
NEWSTR: 


NEWPSUs 


SETPSU: 


END: . 
ELSE IF BYPSU(IV) > BYPSU2(TY) THEN GO TO MEwWPSU 


ELSE IF BYPSUCIV)<BYPSUZ(IV) THEN CALL ERROR(ERRZF ES 


END: 
PROCOERS: 


IF NOBS > © THEN CALL ADDSUMS: 
IF LASTPSU=’C’ THEN DO: 
C&LL PRNTSTR: 
PUT EDIT(’ DUMMY PSU CREATED ’)¢(4): 
END: 
STRATNO=STRATNO+1 + 
ISINDX=ISINDX+1: 
IF SUBSTR(DESVEC,ISINDX,1)=’X’ THEN ISIMDM=i: 
DESERR='0'B: 
LASTPSU=/0': 
DO I¥=1 TO NSTRID: 
BYSTR2(IV)=BYSTR(IW): 
END: 
GO TO SETPSU: 


CALL ADDSUMS: 
IF LASTPSU=’0’ THEN LASTPSU=‘1’: 
ELSE DESERR=/1°R: 


NPSU=NPSU+1: 
DO IV=1 TO NPSUID: 
BYPSUZ(IY)=BYPSUC IV): 


02290 
VEB00 
O2310 
02320 
2300 
O254d0 
02350 
02360 
02370 
02389 
02390 
02400 
o2410 
02420 
02430 
o2440 
02450 
2460 
02470 
02480 
O2490 
O2500 


M2519 
2520 


HR2E20 
O2530 
92540 
02550 
02560 
02570 
02580 
02590 
O2600 
02610 
02620 
02630 
02640 
02650 
92660 
92670 
Y— Ls 

N26B0 
02690 
02700 
92710 
92729 
02730 
02740 
02750 
92760 
02770 
02780 
02790 
92600 
O2B10 
02820 
02930 
02640 
02850 
02860 


@ 
ERIC 


PROCOSS: 


ENDWAR s 
ENDCASE: 


END: 


IF DESFLG“*DESERR THEN DO: 
CALL PRNTSTR: 


CALL P 
END: 


IF DESERR 
PUT EDIT(’ MORE THAN 2 PSUS. THIS ONE SHIPPED’) 


CALL Y 


RNTPSUs 


THEN 


ARK(6 WOT): 


IF MISSING(WGT) THEN DO: 


CAL 
END: 
NOBS=N 


L ERROR(ERRiF,ESTAB): 


O8#S+1: WOTORS=WGTOBS+WGT: 


IF DESERR THEN GQ TO ENDCASE;: 


CALL WV 


ARXC1, INHOLD(1)): 


IF NOMSFLG THEN DO fIV=1 TO NW: 


IF 


MISSING(INHOLD(IV)) THEN C 
NMISSSNMISS+13 
WGTMISS=WGTMISS+NGT: 

GO TO ENDCASE: 


ct 


END: 


END; 


In 


=NREP+1: 


KV=03 
IF \NOMSFLG THEN KV=2"NVi: 


ie 


2ENVIFKYs 


DO IV=1 TO N13; 


IF IV = NVi THEN X=1.0: ELSE X=INHOLDC(IY 
KWY=KVU4+1s 
IF \NOMSFLG THEN DQ: 
IF MISSING(X) THEN GO TO ENDJAR: 
ESTAB(IR, IN) SESTABCIR, IV)+1.03 
ESTABCIR, IVAN) SESTAB(IR, IV4ANGL)+HGT: 
END: 
“MWNTEXNEWNGT: 
ESTAB(IR, KY) =ESTABC IR. KY) +KWTS 
ESTABCIR  KV4+NW? )SESTABCIR KVENYS D+KHRWT 
IF IV * 1 & SSQFLG = ‘’1’°H THEN DO JWeil TO 


K=K+1: 
ESTAB(IR,K)=ESTABCIR, K)+NWTS#INHOLD (JY) 
END: 


END; 


END; 


CALL ADDSUMS: 


PUT SKIP EDI 


NUMPSU=NPSU; 


TCNUMPSU,’ PSUS SPECIFIED’, 
NPSU,’ PSUS FOUND’) (F(19),A. SKIP. FC19),A 


iA) 


\a 
ao» 


y; 


/# t/ 
/* IF QUTPUT OPT IN EFFECT CLOSE THE DS.. x/ 
/* #/ 
IF GUTFLG THEN DO : 

CALL SETDSN(2) : 

CALL FRUF +: 

CALL CLSOUT : 
END ; 
ft #/ 


<) 
ERIC 


02870 
QO2BBY 
O2890 
02700 
U29LO 
G2920 
02930 
OSG40 
92950 
O2960 
O2970 
O2980 
92990 
Q{I0o0o 
O3010 
N2020 
OBC3oO 
IBO4O 
O2050 
Q3960 
YS070 
Q23B030 
3090 
23100 
OO1i6 
03120 
YG1C0 
03140 
Q3150 
03160 
03170 
03180 
03190 
GBSo0a 
93210 
Q32Z20 
M3230 
QO3240 
O3250 
03260 
93270 
Q3280 
3290 
03300 
03310 
03320 
O33S 

03340 
03350 
03360 
03370 
Q3390 
O2390 
03400 
03410 
03420 
023430 
03440 
O34E0 
03460 


/* 
PREPOUT:PROC 


7s 
/% 
/* 
/%* 
/# 
/* 
/* 
/* 
/* 


Ln ee 


/% 
/* 
/* 
DCL 


/* 
/* 
/* 
DCL 


/* 
/* 
/* 
DCL 


/* 
/% 
/%* 
CALe 
CALL 
CALL 
CALL 
CALL 
CALL 


PREPOUT-INTERNAL PROC TO DEFINE OUTPUT 2D 
& ALL OUTPUT VARIABLES.. 
OUTPUT DS WILL CONTAIN THE 


uw 


FOLLOWING: 

21. REPLICATE ID REPLID 

2. ANY BY WARS DEF ON INPUT 
3. REPLICATE ESTS DEF CN INPUT 


WAR DEF STRUCT.. 


1 NAMESTR, 
NTYPE FIXED BIN(1S), 
NPOS FIXED BIN(15), 
NLNG FIXED BIN(1S5), 


NYWYARO FINED BIN(15), 
NNAME CHAR(S). 
NLABEL CHAR( 40). 
NFORM CHAR(&?, 
NIFORM CHAR(®), 


NFL FIXED BIN(15), 
NFD FIXED BINC1S), 
NF FIAED BINCLIS? , 


NM ph fh) pi fite fin N 


NJUST FIXED BIN(iS) : 
SAS LINKAGE.. 


GPNCUT ENTRY, 

GETDS ENTRY(PTR), 

HISTRY ENTRY (PTR). 

NAMEY ENTRY(FIXED BIN(S1),FIXED BIN(31), 
FIXED BIN(15),), 

ONAME2 ENTRY(FIXED BIN(15) FIXED BIN(3i)). 
ONAMES ENTRY(FIKED BIN(1S) FIXED RIN(31), 
FLOAT BIN(S3)), 

ENDNAM ENTRY, 
DLTOUT ENTRY, 
ABUF ENTRY 3; 


OUTYAR STRUCTURE... 


OUTSTR(1) BASED(OUTPTR). 
OUTADDR PTR, 

VARADDR PTR, 

OUTNAME CHAR(B) : 


NNN Fe 


DEFINE OUTPUT DS.. 


SETDSN(2) +: f/* OUTPUT IS CURRENT. 
OPNQUT : /* OPEN QUTPUT Ds.. 
SETDSN(1) 3; 7* INPUT IS CURRENT.. 
FBUF ; /* FREE INPUT BUFFER 
GETDSCIPTR) : /* ADDRESS INPUT.. 
SETDSN(2) : /* OUTPUT IS CURRENT 


x / 
%/ 
#/ 


#/ 
Hf 


a/ 
#/ 


03470 CALL. RISTRYCIPTR) + /* GET DS HISTORY.. w/ 


o34g80 /% ay 
03490 /* DEFINE REPLICATE IDENTIFIER.. Hy 
O3800 /% t/ 
O3510 NNAMES‘REPL_ID’ + 

Ossz20 NLABEL=’REPLICATE IDENTIFIER? : 

02530 NFORM.NIFORM=’ ‘ ¢ 

o3540 NTYPE=1 : 

lc fatat®) NLNG=8 : 

03560 NFL,NFD= 0 + 

03570 CALL ONAMES(NTYPE.OUTYVEC, REPL ID) 3 

o3580 /# % / 
03590 /* DEFINE VARS ON VARLIST.. a / 
03400 /* %/ 
03610 DO IV=1 TO Nui : 

034620 CALL SETDSN(1) ; 

03630 CALL NAMEV(1,IV.NTYPE,NAMESTR) 

03640 NLNG=S : 

034650 CALL SETDSN(2) 

03660 CALL ONAMES(NTYPE,OUTWEC,FLAT_EST(1V)) 

OBE7O END : 

o268e Js #/ 
o2450 /# DEFINE COMPUTED VARS.. #/ 
O2 7OO /* * / 
Q3710 IF OUTVFLG THEN DO IV=i TO NSTAT: 

03720 NNAME=OUTNAMEC( IV) : 

03730 NLABEL=’COMPUTED ESTIMATE’ : 

03740 NLNG=8 + 

O37E0 NTYPE=1 : 

93760 NFL .NFD=0O : 

93770 NFORM, roee ‘os 

03780 CALL ONAMES(NTYPE, OUTVEC .FLAT_LEST(IVENW1)) 
03790 END : 

O3800 /% “/ 
03810 /* TERMINATE DS & OUTPUT WAR DEFINITION.. #/ 
03820 /* DELETE OUTPUT ROUTINES... / 
03830 /%* é 
02840 CALL ENDNAM : 

O8B850 CALL DLTOUT : 

03860 CALL ABUF : 

03870 CALL SETDSN(1) : 

o3880 CALL ABUF : 

03890 END PREPOUT : 

C3900 ADDSUMS: PROC: 

OBE10  f* INTERNAL PROC TO ADD PSU SUMS +/ 

3920 /% TO HALF-SAMPLE SUMS #/ 

03930 IRIZNREP+1% 

03940 WGT=1.0: 

03950 IR INDX=ISINDX: 

03960 DO IR=1 TO NREP: 

03970 IF IR > 2 THEN Da: 

02980 IRINDX=IRINDM+1: 

a3990 ' TESTPSU=SUBSTR(DESVEC, IRINDM.1): 

04000 IF TESTPSU = ’X’ THEN DO: 

o4o10 IRINDX=i1: TESTPSU=SUBSTR(DESVEC,1,1)5 
04020 END; 

04020 IF DESFLG THEN IF LASTPSU=TESTPSU 

04040 THEN PUT EDIT(’O ‘) (A): 

04050 ELSE PUT EDIT(‘’1 ’) (A): 


<) 
ERIC 


O4060 IF LASTPSU=TESTPESU THEN GO TO ENDREP: 


04070 END: 

04080 7* SECOND REP INCLUDES LAST OF EACH PATR 4/ 
N4090 ELSE IF IR = 2 THEN DO: 

04100 iP PESFLG THEN PUT ERIT(LASTESUL’ “SCALA: 
n4i10 WGT=2.08HGT: /# DOUBLE WEIGHT FOR HALF EFAMp 
LES #/ 

o4220 IF LASTPSU = ‘0’ THEN GO TO ENDREF: 

04130 END: 

04140 DO IV=i1 TQ NSTAT: 

04150 ESTABCIR, IN) =ESTABCIR, IV) + 

04160 WGTSESTAB( IRL, iy): 

04170 END: 

04180 ENDREP: END; 

04190 DO IW=s1 TO NSTAT: 

M4200 ESTABCIR1,IV)=0.0: 

04210 END; 


04220 END ADDSUMS: 
04230 PRNTSTR: PROC: 


M4240 /% INTERNAL PROC TO PRINT STRATUN ID af 
O4250 PUT SKIPs 

04260 DO IVsi TO NSTRID: 

04270 SF BYCHARZ( IV) THEN 

04280 PUT EDIT (BYCHSTZ(IW).’ ()(A,A): 

04290 . ELSE PUT EDIT (BYSTREC(IW).% “) (FCG) Ad: 
04300 END: 

04320 END PRNTSTR: 

O4320 

C4530 PRNTPSU: PROC: 

04340 /* INTERNAL PROC TO PRINT PSU !D / 

C4350 DO IV=1 TO NPSUID; 

04360 fF BYCHAR( IW) 

94370 THEN PUT EDITCRYCHPS2(IM),’ “)CA,AD: 
04380 ELSE PUT EDIT(BYPSUZ(IV).* 4) (F(B),A)s 
04390 END: 


C4400 END PRNTPSU: 

04410 END PROCESS ; 

04420 /#MEM&#/ 

04430 /# . 

04440 //LKED.SYSLMOD DD DSN=NCES.MEJQNY.SAGELIB.DISP=SHR 
04450 //LKED.SYSLIB DD 


Q4460 // DD 

04470 // DD 

04480 // DD 

04490 // DD 

04500 // DD DSN=NCES. XEJQNY.SAGELIB,DISP=SHR 
043510 // DD DSN=SYS2.SAS.795.SUBLIB,DISP=SHR 


O4520 // DD DSN=SYS2.SAS.R795.GLIBRARY ,DISP=SHR 
043530 //LKED.SYSIN DD # ‘ 

043540 NAME PROCESS(R) 

O4550 /# 

END OF DATA 


re) 
ERIC 


QO010 
Oon2zY 
OOO3 o 
OnoO40 
Odds 
PAGO 
OOo 7a 
ono0g8a 
ONO 
O10 
OO1L0 
OoOLz0 
O0130 
00g140 
QO150 
00160 
00170 


eoLeo 
OOL9G 
OO20C 
OO210 
OO22 

Q9O2Z30 
N0240 
QOO2Z50 
OO2Z60 
CO270 
og2zeaa 
00290 
QIBOAY 
OO910 
00929 
N0330 
O34 
oOogsa 
QO0360 
00370 
M0380 
90390 
ON400 
00410 
00420 
00430 
00440 
90450 
00460 
00470 
00480 
90490 
00500 
09510 
O0520 
If} ) 

00530 
OoOS40 
OP2OSSa 
GCOS450 
00870 


re) 
ERIC 


//K JOR (ED.AIR), ’D-WYCUNG’ .PRTY=8, 

// MEGCLASS=P ,NOTIFY=XEJQNY ,NMSGLEVEL=(2.4) 

7*ROUTE PRINT R207 

f/NYYSTEP2 EXEC MPLIMCL SYSOUTE=P, 

// PARM.PLI=’NM,ND,A.NEST.X,NC(S) ,MAR(2,72,1).NSYN(C).NOP.QF.C’ 
// PARM.LAED= ‘MAP .MREF ,.LIST.NCAL’, 

// COND.LKEDS(4.LT,PLI) 

//PLI.SYSIN ODD *# 


PRINTIT:sPROC(ESTAB) 3; 
/* tt / 
/* EXTERNALS.. #/ 
f% a/ 


DCL (NY1.NOUT,NWGT,NOHS,NREP,NDEP,NMISS.NSTAT) 
FIXED BIN(31) EXTERNAL, 
(NUMPSU WGTOBS “IGTMISS) FLQAT BIN(S3) EXTERNAL. 
(MISSPTR.,MWGTPTR) PTR EXTERNAL, 7 


(SUMFLG.MEANFLG.STDFLG, SSQFLG.COVFLG.COSFFLO.RESRFLG 

NOMSFLG,DUMFLG) BIT(1) EXTERNAL 
/* #/ 
7% INTERNAL STRUCTURES... it / 
/* xy 
OCL ESTAS(#,*) FLOAT BIN(5S3) CONTROLLED, 

1 NAMESTR, 

2 NTYPE FIXED BIN(19), 

2 NPOS FIXED BIN(15). 

2 NLNG FIXED BIN(15). 

2 NVARO FIXED BIN(15). 

2 NNAME CHAR(8?. 

2 NLABEL CHAR(40), 

2 NFORM CHAR(8), 

2 NIFORM | CHAR(8), 

2 NFL FIXED BIN(15). 

2 NFD FIXED BIN(15), 

2 NF FIXED BIN(15). 

2 NJUST £FIXED BIN(15) : 
DCL SYSPRINT FILE PRINT: 
/* %/, 
/* DECLARE LINKAGE TO SAS... 2 / 
/* #/ 


DCL NAMEY ENTRY(FIXED BIN(31).FIXED BIN(31), 
FIXED BIN(15),). 
BEST ENTRYC(FLOAT BIN(S3).FIXED BIN(31), 
FLOAT BIN(S3)?. 
IOPT ENTRY(FIXED BIN(31)) 
RETURNS(FIKED BIN(31)), 
LINSIZ ENTRY(FIXED BIN(31)), 


STITLE ENTRY(FIXED BIN(31),FIXED BIN(31)), 

ERROR ENTRYC(CFIXED BIN(G1).(0#.*#) FLOAT BIN(S3)) : 
/* #/ 
/*  LQCALS.. ‘ #/ 
/* #/ 


Dou (2, IR. IV. IU.KLNV.L.CENTER,LNSZ. TER. KU KY. TENE NPD, 


FIXED BIN(31). 
(ADDR. SUBSTR.SQRT) BUILTIN, 
(EST.SE.LB.UB,KNR.XM,KN.XS) FLOAT BIN(S2), 
FCHAR CHAR(19) ALIGNED, 


DEEST FLOAT BIN(S3) BASED(ADDR(FCHAR? ) . 


an 


oo0seo 
00590 
NO600 
20610 
oné62n 
o0e30 
00640 
OG6E0 
00650 
ae 
20670 
0G680 
90690 
00700 
OC710 
£19720 
co720 
) 
00740 
00750 
00740 
00770 
00780 
QATIDO 
mogoo 
ABEL / 
0810 
RA, 
fOB2Z0 
NOBIO 
OOB8Ad 
aogso 
O0860 
Ne870 
ER’. 


-. 00880 
' - 00890 


ra) 
ERIC 


on900 
>, A, 
00910 
00920 
00930 
00940 
09950 


00960 
ER’, 

90970 
00980 
90990 
>, A, 
01000 
01010 
92020 
O1039 
O1040 
01050 
01060 
01070 


~~ 


f+ 


/% 
{% 


/% 
Dc 
DC 


STATS(10) FLOAT BIN(S2): 


4/ 
ERROR MESSAGES... 4/ 
: 4/ 


L SUSTIT CHAR(32Z) STATIC INIT( ‘SUMMARY INFORMATION’): 
L STATLEL CHAR(10) STATIC: 


DCL ERR1 CHAR(80) INIT( 


/* 
/% 
/% 
ON 


EN 


/* 
CA 
IF 


’ ERROR: LINE SIZE MUST BE > 132 IN PRESENT VERSION 
ERRIF FIXED BIN(31) BASED(ADDR(ERR1)?}: 
DEFINE ENDPAGE CONDITION.. ‘ x/ 


ENDPAGE(SYSPRINT) BEGIN : 
CALL STITLE(O,L) ¢ 
PUT EDIT( ‘PROC BRRVAR — BRR SAMPLING ERROR ESTIMATES’ 


(SKIP (1) .COL(CENTER) .A) : 
PUT SKIP EDIT(SUBTIT) (COL(CENTER) ,A): 
PUT SKIP (2) : 
IF SUBSTR(SUBTIT.1,1)» = ‘U’ THEN DO: 
PUT EDIT( ‘95% CONFIDENCE BOUNDS’, ‘STATISTIC’. 
‘VARTABLE’, ‘WTD. N’%, 
‘ESTIMATE’, ‘STD. ERROR’, ‘LOWER’. ‘UPPER’. ¢L 


(COL(468).A,SKIP,A,COL(12) ,A,COL(26) ,A.COL(3S?. 


COL(S0) .A,COL(67) .A.COL(80) .A, COL (87) Ad: 
END: 
SLSE IF SUBSTR(SUPTIT,.1,1) = ’C* THEN DO: . 
PUT EDIT(’*FIRST’,.. “SECOND’. SUBSTR(SUBTIT,£.11).. 
‘95% CONFIDENCE BOUNDS’. ‘WARIABLE’. ‘WARIABLE 


‘WTD. N’. ‘ESTIMATE’. ‘STD. ERRGR’ “LOWER °, ‘UPP 


“PROB C>0’) 
(A. COL(12), A, COL(22), A, COL(63), 4. SKIP, 
A, COL(12). A, COL(26). A, COL(29), A, COL¢(S1 


COL(64). A. COL(77), A.COL(87).A): 
END: 
ELSE IF SUBSTR(SUBTIT,1,1) = ‘R’ THEN DO: 
PUT EDIT(‘’FIRST-. ‘SECOND’, SUBSTR(SUBTIT,1,11), 
“'95% CONFIDENCE BOUNDS’, ‘VARIABLE’, ‘WARTABLE 


“WTD. N’,“ESTIMATE’.’STD. ERROR’. ‘LOWER %. ‘“UPP 


‘PROB C>0’) 
(A, COL(12), A, COL(22), A, COL(68), A. SKIP. 
A, COL(12), A, COL(26), A, COL(39). A, COL(S1 


COL(64), A. COL(77), A, COL(87)_A): 
SND: 
PUT SKIP: 
D; 


SET LINESIZE AND CENTER +f 
LL LINSIZ(LNS2Z):3 : 
LNSZ < 80 THEN CALL ERROR(ERRIF. ESTAR): 


01080 CENTER=(LNSZ-35)/2: 


01090 IF IOPT(SO) = 0 THEN CENTER=1: 

01100 

oi1190 : /* PRINT SUMMARY REPORT af 
O1120 SIGNAL ENDPAGE(SYSPRINT): 

01130 PUT SKIP SDITCNOQBS,’ OBSERVATIONS FESD’) 
0114060 (F(12),A): 

01150 PuT SKIP EDIT(WGTOBS. ’ WEIGHTED GESENVATIGNS READ’) 
O1LLS¢ (F(12.1), Ads 
01170 PUT SHIP: 

02180 IF NOMSFLG THEN DO: 

QO1190 WGOTCBS=WGTOBS-ESTABCIINUI): 

02200 1 PUT SKIP EDIT(WGTOBS. 

O1219 ‘ WEIGHTED OBS DELETED FOR MISSING VALUES’) (F (12, 
1),Ads 

01220 END: 

01230 ELSE PUT SKIP “EDIT(’ MISSING 4ALUES DELETED SEPSRATELY 
O1240 “FOR EACH VYAR‘’ICALAD: 

O1250 PUT SKIP: 

22260 PUT SKEP(2) EDIT(NREP,’ REPLICATES IN’ LNUMPSU, 
OLE 70 “PSU PAIR DESIGN’) (F(12),%(1),.0,F (4) 4615.0): 
©1280 

01290 /# COMPUTE AND PRINT Sums %/ 
91300 SUBTIT = ‘’UNIVARIATE STATISTICS’: 

01310 NYVENVYG—t:  XNR=NREP-1: 

01220 IF NOMSFLG THEN KVi=0; 

012336 ELSE KY¥L=2HNY1: 

O1G46 IF SUMFLG THEN DO: 

91350 SIGNAL ENDPAGE(SYESPRINT?: 

71360 STATLBL=’SUM‘: 

Q1370 XNS=ESTAB(CIINYU1)s 

01380 KV=SKY 1s 

01390 DO IV=1 TO Nis 

01400 KYSKW+ls 

014106 IF \NOMSFLS THEN KN=ESTAB(1,IV4NU1): 
01420 CALL PRINTLNCKW, IVLXN) s 

01430 PUT EDITC(NLABEL) (COL(74).4): 

01440 END: 

01450 END: 

01460 

O1470 /*# COMPUTE MEANS Ht / 
01480 DO IR=1 TO NREP: 

01490 KVEKVU1s 

01500 XNFESTABCIR, NVI): 

01510 DO IV=1 TO NY: 

01520 KUY=KY+4 2 

01530 IF \NOMSFLG THEN KN=ESTABCIR, IU4NU1): 
01540 ESTABCIR. KV) =ESTABCIR.KY)/XN: 

O1550 END: 

01560 END; 

01570 

015280 IF MEANFLG THEN DO: 

01599 STATLBL=‘MEAN’: 

01600 KN=ESTAB(1LNYI): 

01610 KYV=KW 13 

O1620 DG IV=1 TO NV: 

01630 KYSKU+i = 

01640 IF \NOMSFLG THEN XN=ESTAB(1,IV4ENU1) = 


01650 CALL PRINTLENC(KY, IV, XN): 


° 
ERIC 


91660 PUT EDIT(NLAPEL) (COL(87) 43s 


01670 END: 

1680 END: 

1690 

01790 /* COMPUTE S. D.S #/ 
O1710 Po ik=1 TO NREP: 

01720 MNSESTABCIR NVI): 

01730 KVEKW1s 

01740 DO IV=1 TO NV; 

o1L750 KVEKW4is 

01740 IF \NOMSFLG THEN XN=ESTAR(1,IV+NV1): 
01770 SE=ESTABCIR.KV+NV1)/XN-ESTASLIR, KY) #2: 
o17BC IF SE > 0 THEN SE=SQRT(SE): 

01790 ESTABRCIR, KV+NV1)=SE: 

o1200 END: 

01816 END: 

01820 : 

Y1B30 IF STDFLG THEN DO; 

01840 STATLBL=’STD’: 

o1es0 . XN=ESTAB(1,NV1): 

01860 KWEKW1: 

01870 DO IV=1 TO NW: 

01990 KV=EKV41s 

o18e0 . IF \NOMSFLG THEN KNSESTAB(1, TN4NY1) + 
91900 . CALL PRINTLNEKVENYL, IY SXN): 

01910 PUT EDIT(NLABEL)(COL(74).A): 

01920 END; 

01930 END: 

‘01940 

01950 /* COMPUTE COVARIANCES AND CORRELATIONS #/ 
2960 IF SSQFLG & NY > 1 

01970 THEN DO IR=1 TO NREP: 

01980 K=2eNV1s 

01990 DO I¥V=2 TO NY: 

2900 XM=ESTABCIR, IV): XNSESTGB(IR, NVI): 
C2ZOLO ‘ DO JV=1 TO IV-1: 

02020 K=K+13 

02030 | ESTABCIR,K)=(ESTAB(CIR.K)/XN-ESTABCIR, JU) #XM) + 
02040 END: 

92050 END: 

02060 END; 

02070 

02080 | IF CORRFLG THEN DO: 

02090 - DO IR=1 TO NREP: 

02100 K=2#NV13 

02110 DO Iv=2 TO NW: 

02120 XS=ESTABCIR, IV+NV1); 

02130 DO JV=1 TO I-21: 

02140 K=K+1: 

02150 ESTAB(IR.K)=ESTABCIR,K) 
02160 /(XS#ESTABCIR, JV+NV1)): 
02170 END: 

02180 END: 

o2190 END: 

02200 END: 

C2210 . 

02220 IF COVFLG*CORRFLG THEN DO: 

02230 IF CORRFLG THEN 

02240 SUBTIT=/CORRELATION SSTIMATES’: 


02250 ELSE SUBTIT=’COVARIANCE ESTIMATES’: 


ERIC 


C2260 ANFESTABCILNYI): 


92274 SIGNAL ENDPAGE(SYSPRINT): 

92280 K=24nis 

02290 , DO IV=2 TO NY: 

02300 CALL NAMEV(1, IVINTYPE NAMESTAR): 

M2319 STATLBL=NNAME: 

022320 DO JWV=1 TO IVv-ts 

02330 K=K+1: 

O2540 CALL PRINTLN(K, JY XN): 

O2350 END: 

02560 END; 

02370 END: 

02380 

02390 DCL CPTR.PTR, Ke FIXED BIN(31), 

O2400 C(1) FLOAT BIN(S3) BASED(CPTR): 

M2410 If REGRFEFLG THEN DO: 

02420 SUBTIT=’REGRESSION COEFFICIENTS’: 

Y2450 SISNAL EWDPAGE(SYSPRINT): 

O2440 IRI=aNREP+1: 

E2450 CPTR=ADER(ESTABCIR1I.1)): 

M2460 NP=NVYI-NDEP-1: 

02470 NP1=NP+15: 

M2480 DO If=1 TO NDEP: 

92490 KYSNP+IT: 

O2500 DO IR=1 TO NREP: 

92510 K=24NU1s 

92520 K2=0: 

o2530 DO I¥=1i TO NP: 

e2540 IF IV > 1 THEN 

O2550 DO JV=1 TO Tvei: 

Q2Z5E0 K=K+15 K2=K2+1: 

C2570 C(K2)=ESTABCIR.K): 

O2580 END: 

92590 K2=K2+1:3 

C2600 IF CORRFLG THEN C(K2)=1.0: 

02510 ELSE C(K2)=ESTAB(IR, IV4tNW1) #2: 
_ 02620 END: — 

02630 K=(KYU-1)#(KU—2) /2+22NU 1s 

92640 DO IV=1 TO NPs 

02650 K=K+i: K2=K2+1: 

02660 COCK2Z)=ESTABCIR,K): 

024670 END: 

02680 K2=K2+13 

02690 IF CORRFLG THEN XS=1.0: 

02700 ELSE XS=ESTABCIR, LI+NP+NV1) #2: 

02710 C(K2)=XS:5 

M2720 * CALL REGRLTS°,NP1,I1ER): 

02730 , K=O: 

027409 K2=NP#(NP+1)/2: 

02750 DO IV=1 TO NP: 

02760 K=K+13 K2=K2+1: 

027706 ESTABC IR. KJ =C(K2) 2° 

02780 END: 

92790 K=K+12 K2=K2+15 

02800 ESTABCIR.K)=(XS-C(K2))/XS: 

02810 END: 

02820 CALL NAMEV (1. KYNTYPE,NAMESTR): 

O2830 STATLBL=NNAMEs 

02840 DO IV=i TO NP: 

02850 CALL PRINTLNC IVY, IV.XN): 


° 
ERIC 


22ZB60 
O2B70 
o2esa 
H2aS0 
92900 
g2Fid 
YZ920 
02930 
02940 
2950 
N2946C 
297% 
02980 
C2990 
' aged 
O3010 
02020 
03030 
OSO40 
OBOE0 
O3060 
OOo7a 
82080 
V3090 
03100 
03110 
O3120 
O3130 
OBi40 
O3150 
L) 
03140 
“03470 
O2286 
03190 
OBZ90 
03219 
Ria / 
03220 
03230 
03240 
03250 
03260 
03270 
032990 
03290 
033500 
03310 
03220 
03330 
03340 
03350 
03369 
03370 
033890 
03290 
02400 
03410 
C3420 
67 430 


ERIC 


END: 
STATLEL=/MULT RSQ’: 
CALL FRINTLNONPL KY MND: 
END: 
END: 
/* PRINT ESTIMATES FOR EACH RALF SAMPLE af 
/* s/ 


IF BDUMFLG TREN DO: 
SUBTIT=‘REPLICATE STATISTICS’: 
SIGNAL ENDPAGE(SYSPRINT): 
DO IV=1 TO NSTAT: 
PUT SKIP: 
DO IR=i TO WREP: 
PUT ENITCESTABCIR. IW) (CF (8,2)): 
END: 
END: 
END: 


PRINTLN:s PROCCK.IWV XN); 
/* SUBROUTINE TO SET UP AND PRINT & SINGLE LINE #/ 


/*® K POINTS TO THE ESTIMATE IN ESTAS af 
DCL (K,IW) FIXED BIN(31): 
DCL XN FLOAT BIN(S3): 


CALL NAMEY(L.IY,NTYPE.NAMESTR) 

EST=ESTAB(1,K); 

CALL SAMPERR(ESTAB. K, XNR, SE): 

LESEST-1.96#SE;3 UB=EST+1.74#SE; 

PUT SKIP EDIT(STATLBL, NNAME. XN. EST. SE. LE. UB. NLABE 


fA, COL(12). A, COL(20), F(19,1).F(19.2).F (19.9). 
PLIG,3). F(ia.3), Ads 
END PRINTLN: 


SAMPERR: PROC(ESTAB, II, XNR, SE); 
/# INTERNAL SUBROUTINE TO COMPUTE THE SAMPLING STD ERRO 


DCL ESTAB(#, *) FLOAT BIN(S3) CONTROLLED; 


DCL (SE, XNR) FLOAT BIN(S3), 
(IR, II) FIXED BIN(31), 
NREP FIXED BIN(31) EXTERNAL: 
SE=0.0: 


EST=ESTAB(1,II); 
DO IR=2 TO NREP; 
SE=SEt(ESTABCIR, II)-EST)##2: 
END: : 
IF SE > O THEN SE=SQRT(SE/XNR) : 
END SAMPERR; 
REGRLT: PROC(C,NV IER): 


/%* REPLACES A LOWER TRIANGULAR COY MTX %/ 
/%* WITH A MTX OF REGR COES WITH ERROR x/ 
/* WARS ON THE MAIN DIAG. (EACH VAR IS Hf 
/* PREDICTED BY THE PRECEDING ONES) %/ 
/%* IER IS SET TO THE NUMEER #/ 
DCL C1) FLOAT BIN(S3), 


{B,D.E) FLOAT BIN(S3) STATIC, 
(NY IER) FIXED BIN(31), 


yd 


03440 CN,KRLI,KP,KR,AC KE, JJ.0) FIXED BIN( 31) STATIC: 


03450 

O3460 TER=0; 

03470 N=NY—j: 

03480 KP=0: 

03490 DO K=1 TO WN: 

O3500 KR=KP?: KP=KP+Ks 

o3s5i0 D=C(KP): 

03520 TF DPD < 1.E-20 THEN Do: 
03530 TERSi: 

O3540 GO TO ENDREG: 

C2550 END: 

O3560 DO I=1 TQ NW: 

03576 IF I < K THEN DO: 
023580 KR=KR+12 KCSKP+Ks JJ=K+15 
03590 END: 

036c0 ELSE IF I » K THEN DQ: 
03410 KR=KR+I-i: KC=KR: dJ2ts: 
03620 END: 

o2e3o ELSE DO: 

02640 KRSKR+1: 

03650 GO TO NEXTI: 

03660 END: 

02679 S=C(KR)/D: 

024580 DO .J=JJ TO NW: 

03490 KESKC+I-K: 

03700 C(KE)SCCKE)-S#O(KC): 
O3710 KCSKC+J3 

2720 END: 

03730 IF I >» K THEN C(KR)=B: 
62740 NEXTI: END: 7 ; 
O3750 END: 

O3760 ENDREG: 


03779 END REGRLT: 

03780 END PRINTIT: 

03790 /* MEMB #/ 

OGB00 /* 

Q3E10 //LKED.SYSLMOD DD DSN=NCES.XEJQNY.SAGELIB.DISP=OLD 
OGB20 //LKED.SYSLIB DD DSN=NCES.XEJQNY.SAGELIB, DISP=SHR 


03830 // DD DSN=SYS2.SAS.V796.SUBLIB.DISP=SHR 
03840 // DD DSN=SYS2.PLIXBASE.DISP=SHR 
03850 // DD DSN=SYS2.PLIKLINK, DISP=SHR 


03260 //LKED.SYSIN DD # 
03870 NAME PRINTITCR) 
NSBB0 /# 

END OF BATA 


@ an 
ERIC 


@ 


ERIC 


madg10 
M3020 
2C030 
90040 
Qo0080 
00060 
OOa7a 
00080 
HoOO90 
00100 
ogo110 
oeo1z0 
00130 
OOL40 
OO1Sa 
N01E0 
00170 
oo1ea 
OL90O 
0200 
O0210 
9220 
o230 
vO2Z40 
90250 
ug269 
Q0270 
9o280 
C0290 
00300 
90310 
O0320 
OO330 
00340 
OO350 
02360 
00370 
00380 
vO3IO 
00490 
00410 
00420 
90430 
00440 
00450 
eo460 
00470 
o0480 
90490 
Qe300 
90510 
00520 
OOS520 
00340 
003550 
oas60 
00570 
o058¢ 
00590 


//& 
// 


/*ROUTE PRINT 
7/ASM EY 


{EC 


PARSER. 


JOB (ED,ATR). 7D-WYOUNG’ .PRTY=8. 


MSGCLASS=P ,NOTIFY=XEJOQNY .MSGLEVEL=C(1 1) 


ROD 
MASMYCL , CREGSS500K 


//ASM.SYSPRINT DD SYSOUT=P 


//ASM.SYSLTB DD DSN=8) 


//ASM.SYSIN DD # 


PROC 
OPTS 
OPTS 
Q°TS 
OPTS 
OPTS 
OPTS 
OPTS 
OPTS 
OPTS 
PARMS 
PARMS 
PARMS 
LIsTs 
LISTS 
LIsTs 
LISTS 
LISTS 
LIsTs 
LISTS 
DATASETS 
DATASETS 
DATASETS 
DATASETS 
DEF SULT 
DEFAULT 


INIT 


INIT 
PSUEDO 


BRRVARD 
COMPLINK 
BRRVAR 
TERMER 


TERMSUB 


TERMSUB 


PRINT 
SASPROC 
SASLIST 
SASLIST 
SASLIST 
SASLIST 
SASLIST 
SASLIST 
SASLIST 
SASLIST 
SASLIST 
SSSL Ist 
SASLIST 
SASLIST 
SASLIST 
SASLIST 
SASLIST 
SASLIST 
SASLIST 
SASLIST 
SASLIST 
SASLIST 
SASLIST 
SASLIST 
SASLIST 
SASDFLT 
SASDFLT 
BR RT 
ENTER 
SASLOAD 
BALR 
STM 

MVC 

BAL 


SASDLT 
LEAVE 


NOGEN 


NAME>=BRRVUAR. LOADMOD=ERRVAR?, 
MEAN.1,MEANS .1 


STD.2 - 
CoV ,3 


(‘o? 


[woe 


COR.4,CORR.4 
REG,S,REGR & 
SUM,4.SUMS 16 


DES.7,DESIGN.7 


DUMF .8 
NOMISS 


MUMPSU.1 NPSU,1  MODE=NUMERIC 
NRES ,2.NRESS 2. MODE=NUMERIC 
NDEP .3,.NDEPS,2,MODE=NUMERIC 


~9,NM.9 


SAS.4796.MACLIB _ DISPESHE 


DATA=ASGIVEN 


VARISELES.1,VAR,1,MODE=NUMERIC 
CLASSES ,2.CLASS,2_.MODE=ALFPHANUM 


PSUID.2.PSU.3,.MODE=ALPHANLUM 
ADJBY,4.MOQDESALPHANUM 


STRATID.S.STR.5,MODE=ALPHANUM 


WEIGHT.6,WGT,6.MODE=NUMERIC 


COMPUTE ,7,.COMP.7,MODE=NUMER IC 


DATA, 1.MODE=DSIN 
OUT.2,MODE=DSOUT 
OUTEST,3,MODE=DSOUT 
ADJDATA,4,MODE=DSIN 
NO=1.NAME=_LAST_,Ma@n==DSIN 
NO=2,NAME=_DATA_,MODE=DSOUT 


SASART 
14,15 


O,1,COMPLINK 


@NAME.=CLE’_OBS_ ’ 


RT .PSUEDO 


@MAJNAME.C’ * 


CSW, @SWPSE+@SNDEF 


BMODE ,@A+@S 
15.¢ADDLOC 


iS 


2F 


15.0, COMPL INK 


1,1 
15 


15,0,COMPL INK 


O0600 STMTREC ENTER 


agé10 _ Comp OUTYAR _OUTVAR 
0620 COMP OUTVAERS. QUTVAR 
00630 _ Comp OUTY .OUTYAR 
99640 compe COMPYAR.CMPVAR 
OOES0 COMP COMPY_CMPVAR 
00660 comp COMP  CMPYVAR 
OCcs7e LM 15.0.COMPL INK 
Océse SR Let 

00690 BALR 14.15 

00700 STMTREC LEAVE 

Go710 OUTVAR La 0.16. 

C0720 FAL RT. STMTST 
00730 WORD : 

00749 SR 3.5 

09750 OUTL comp C’:!,OUTE 
00760 LIST OUT! .OUTL 
Oo77E NAME 

QO789 DATA DATAEHN1 N= 
OOFEO LA 5.1(5) 

aral=t81 8) B OUT! 

NOB1iG QUTES ST 5.SUBLIST 
QOLZO SAL RT. STMTEND 
oog830 'B - STMTRECKX 

GOO840 CMPYAR L 6.8DSADDR 
OOBSO USING DSTEL.6 

309860 WORD . 

OOB70 CMPL COMF C’:’.CMPE 
20880 LIST COMPI.CMPL 
20890 LA oO, DSV4R 

rarer-rele) ST 0,@LIST 

OOF1O : D0 SVARF IND 

angz0 LH 0, @VARNO 

00930 LTR 0,0 

00940 Bz COMP I- 

00950 LPADD DSVARLST+i6 
OFS MYC LNAME(12), @NAME 
D097G MYC 22(2,2).8VARNO 
00980 MUT @MODE.84+eS 
DOFEQ MYI @SW,@SWDEF 
01000 MYT @MAJNAME.C’ ’ 
o1010 DO ZADDLOC 

1026 B ; COMPI 

O1030 CMPE - fi STMTRECX 
1040 SASEND INIT=INIT. TTMTREC=STMTREC.TERMERR=ETESMER, TERM=TERMSU! 
01050 LTORG 

01060 END 


01070 //LKED.SYSLMOD DD DSN=NCES.XEJQNY.SAGELIB.DISP=SHR 
O1080 //LKED.SYSPRINT DD SYSOUT=P 

O1090 //LKED.SYSIN DD * 

01100 SETSSI ACOD0001 

O1110 NAME BRRYAR(R) 

END OF DATA 


r) 
ERIC 


